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SUMMARY 


A  mathematical  theory  is  developed  for  the  low-speed  aero¬ 
dynamics  of  rectangular  wings  with  side-edge  separation.  This 
theory  is  then  extended  to  wings  with  sweptback  leading  edges 
(and  straight  trailing  edges)  by  representing  the  actual  wing 
as  a  system  of  elementary  rectangular  wings  of  varying  aspect 
ratio.  Thus,  in  the  limit,  steady  separation  along  the  entire 
leading  edge  is  approximated,  and  the  theory  leads  to  an  iterative 
computational  procedure  for  calculating  the  aerodynamic  charac¬ 
teristics  of  sweptback  wings  with  leading-edge  separation. 

The  assumed  vortex  system  of  the  elementary  rectangular  wing 
consists  of  a  lifting  line  with  its  associated  trailing  sheet 
(lying  in  the  plane  of  the  wing)  and  a  separated  vortex  sheet  from 
each  side  edge  with  an  associated  bound  vortex  system  in  the  wing. 
The  Kutta  condition  is  satisfied  along  the  side  edges,  and  the 
boundary  condition  of  no  flow  through  the  wing  is  satisfied  along 
a  selected  control  line  for  all  spanwise  points  by  means  of  Fourier 
analysis.  The  resulting  theory  includes  the  classical  lifting¬ 
line  theory  of  Prandt]  and  provides  a  means  of  calculating  the 
span  load  distribution  as  well  as  the  strength  of  the  shed  vortices 
and  the  resulting  downwash  throughout  the  flow  field.  The  shedding 
angles  of  the  separated  vortices  are  found  by  an  iterative  tech¬ 
nique  in  which  the  vortices  are  required  to  be  shed  at  the  angle 
of  the  local  flow  immediately  outside  the  side  edges. 

Calculations  are  presented  which  demonstrate  convergence  of 
the  method  for  both  rectangular  and  delta  wings,  both  with  and 
without  iteration  on  the  shedding  angles  of  the  separated  vortices. 
Comparisons  with  experiment  are  also  presented  for  aspect  ratios 
up  to  five,  and  it  is  found  that  for  rectangular  wings  the  cal¬ 
culated  normal  forces  using  one  wing  element  are  in  all  cases 
within  10  percent  of  experiment  but  the  calculated  shedding  angles 
are  unrealistically  high.  For  delta  wings,  the  high  shedding 
angles  evidently  cause  high  predicted  normal  forces.  However, 
the  present  theory  gives  very  good  agreement  for  normal  force  and 


iii 


center  of  pressure  for  delta  wings  of  aspect  ratios  from  one  to 
four,  provided  that  an  appropriate  shedding  angle  is  assumed  and 
held  fixed.  In  this  event,  the  present  theory  gives  substantial 
improvement  over  linear  theory  for  the  prediction  of  normal  force 
and  center  of  pressure,  particularly  at  the  lower  aspect  ratios. 
The  required  shedding  angle  is  found  to  depend  only  on  aspect 
ratio  and  decreases  nonlinearly  with  increasing  aspect  ratio. 
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A  THEORY  FOR  THE  LOW-SPEED  AERODYNAMICS  OF  STRAIGHT 
AND  SWEPT  WINGS  WITH  FLOW  SEPARATION^ 

1 .  INTRODUCTION 

When  a  thin  rectangular  wing  is  set  at  a  small  angle  of  attack 
in  a  uniform  incompressible  stream,  the  flow  remains  attached  to 
the  wing  surface  so  that  all  of  the  shed  vorticity  lies  in  a  sheet 
which  contains  the  wing  planform  and  proceeds  to  roll  up  gradually 
after  leaving  the  wing  trailing  edge  (see  sketch) . 


The  aerodynamic  behavior  of  such  wings  at  low  angles  of  attack  is 
well  predicted  by  the  classical  lifting  line  theory  developed  by 
Prandtl.  As  the  angle  of  attack  is  increased,  the  side  edges 
become  more  oblique  to  the  free  stream,  the  fluid  can  no  longer 
negotiate  the  180°  turn  at  the  side  edges  (from  lower  to  upper 
surface),  and  side-edge  separation  occurs,  giving  rise  to  two 
additional  vortex  sheets  which,  at  higher  angles  of  attack,  begin 
to  roll  up  even  ahead  of  the  wing  trailing  edge  (see  sketch) . 
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Finally,  at  still  higher  angles  of  attack,  the  flow  separates  from 
the  leading  edge,  giving  rise  to  vortices  whose  axes  are  essentially 
normal  to  the  free-stream  direction.  For  this  reason,  the  leading- 
edge  vortices  must  pass  downstream  with  the  general  flow,  producing 
an  unsteady  flow  with  vortex  shedding  in  the  manner  of  flow  past 
a  bluff  body  (see  sketch) . 


The  angles  of  attack  at  which  the  foregoing  phenomena  occur 
and  the  rates  of  rolling  up  and  shedding  depend  upon-  the  aspect 
ratio  of  the  wing,  the  sharpness  of  its  various  edges,  and 
upon  the  Reynolds  number  of  the  free  stream.  In  the  early 
years  of  airplane  flight,  wings  of  practical  design  were -of  rather 
high  aspect  ratio  with  well-rounded  leading  and  side  edges.  Con¬ 
sequently,  separation  from  the  side  edges  was  negligible,  and  the 
classical  lifting  line  theory  was  applicable  over  practically  the 
entire  usable  angle-of-attack  range  up  to  the  onset  o^’  the  familiar 
stall,  corresponding  roughly  to  the  leading-edge  reparation  dis¬ 
cussed  above. 

More  recently,  however,  both  wing  thickness  and  aspect  ratio 
have  diminished  (to  gain  high-speed  performance)  so  that  the  inter¬ 
mediate  regime  of  side-edge  separation  has  become  increasingly 
important.  Furthermore,  for  wings  with  sweptback  leading  edges 
(e.g.,  delta  wings)  the  intermediate  regime  of  steady-state  separ¬ 
ated  flow  is  typical  of  the  landing  condition,  since  such  wings  do 
not  "stall"  in  the  usual  sense.  It  is  in  this  connection  that  the 
present  investigation  has  been  undertaken. 


The  analytical  work  of  the  existing  literature  will  first  be 
reviewed  as  it  bears  upon  the  physical  problem  discussed  here,  and 
then  a  new  theory  will  be  developed  for  the  analytical  treatment 
of  rectangular  wings  with  side-edge  flow  separation.  This  theory 
will  then  be  extended  to  the  case  of  sweptback  wings  with  leading- 
edge  separation,  using  the  rectangular  wing  as  an  elementary  build¬ 
ing  block.  Convergence  of  the  method  will  be  demonstrated,  and 
numerical  calculations  of  normal  force  and  center  of  pressure  for 
both  rectangular  and  delta  wing  planforms  will  be  presented. 
Finally,  comparisons  will  be  made  with  available  experimental  data, 
including  some  span  load  distributions  and  vortex  shedding  angles 
for  delta  wings,  and  the  limitations  of  the  theory  will  be  dis¬ 
cussed  in  some  detail. 

2.  LITERATURE  SURVEY  -  STATE  OF  THE  ART 

The  fact  that  there  exists  an  apparent  fundamental  difference 
between  the  aerodynamic  behavior  of  rectangular  and  triangular 
wings  at  high  angles  of  attack  has  been  recognized  for  many  years. 
It  is  therefore  natural  that  attempts  at  analytical  treatments  of 
high-angle-of-attack  phenomena  have  generally  fallen  into  one  or 
the  other  of  these  two  categories ,  and  only  one  serious  attempt  has 
been  made  to  formulate  a  general  theory  which  would  handle  wings 
of  arbitrary  planform  exhibiting  steady  flow  separation.  Further¬ 
more,  the  inherent  unsteady  character  of  the  flow  at  angles  of 
attack  beyond  the  stall  and  the  strong  dependence  on  viscous  phe¬ 
nomena  have  inhibited  analytical  progress  on  the  problem  of  stalled 
rectangular  wings.  The  first  major  contribution  along  these  lines 
was  made  in  1939  by  W.  Bollay,  vdio  evidently  recognized  that,  for 
rectangular  wings  of  low  aspect  ratio,  separation  occurs  mainly 
along  the  side  edges  and  that  this  type  of  separation  produces  a 
steady  flow. 
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2 . 1  Rectangular  Wings 

2.1.1  Small  aspect  ratio 

The  theory  of  W.  Bollay  (Ref.  1)  treats  rectangular  wings  of 
small  aspect  ratio  under  the  assumption  that  all  shed  vortices  lie 
in  two  planes  normal  to  the  wing  surface  and  containing  the  side 
edges.  Thus,  the  wing  is  represented  as  a  continuous  distribution 
of  horseshoe  vortices  lying  at  some  angle  Q  to  the  wing  surface, 
as  shown  in  the  following  sketch: 


The  angle  0  is  assumed  constant  and  the  vortex  strengths  y  are 
assumed  to  be  constant  across  the  span  (uniform  span  loading)  but 
continuously  varying  in  the  chordwise  direction.  The  chord  load 
distribution  is  taken  to  be  that  of  a  flat  plate  of  infinite  aspect 
ratio  with  a  singularity  at  the  leading  edge;  that  is. 
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where  y  is  the  circulation  per  unit  length  and  c  is  the  chord. 
The  boundary  condition  of  no  flow  through  the  plate  is  satisfied 
"in  the  mean"  along  the  centerline  of  the  plate.  Thus, 
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where  w  is  the  total  velocity  in  the  z-direction  due  to  the 
bound  and  trailing  vortices.  This  procedure  leads  to  a  rather 
complicated  integral  equation  which  is  to  be  solved  for  the  con¬ 
stant  Yq.  Since  the  solution  depends  upon  the  shedding  angle 
6,  Bollay  determines  this  angle  from  the  induced  velocities  by 
imposing  the  condition  that  the  free  vortices  must  follow  the 
streamlines .  This  condition  is  imposed  only  in  the  vertical  plane 
(which  means  the  lateral  velocity  components  are  neglected)  and 
only  at  the  surface  of  the  plate  (which  assumes  that  the  vortices 
in  the  immediate  vicinity  of  the  plate  are  dominant  in  determining 
the  downwash  there).  Thus,  the  shedding  angle  is  interpreted  as 
the  angle  at  which  the  vortices  leave  the  plate  initially. 

Two  limiting  cases  are  considered  by  Bollay  in  Reference  1. 

For  the  case  of  infinite  aspect  ratio  (which  actually  violates  the 
assumptions  of  the  analysis),  the  normal  force  is  found  to  be  inde¬ 
pendent  of  the  shedding  angle  9  and  is  given  by 

Cjj  ■  2Tr  sin  a  cos  a  (3) 

which,  surprisingly  enough,  agrees  with  the  Prandtl  wing  theory. 

In  the  second  limiting  case  of  vanishing  aspect  ratio,  the  shedding 
angle  6  is  found  to  be  just  half  the  angle  of  attack  (for  small 
a) ,  and  the  resulting  expression  for  the  normal  force  coefficient 
is 

C„  -  2  sin^  a  (4) 


i 


which  is  identical  with  the  result  of  Newtonian  impact  theory. 

The  same  result  was  also  obtained  by  Betz  (Ref.  2)  in  1935  by 
introducing  the  "cross-flow  drag"  concept  and  using  the  experi¬ 
mental  drag  coefficient  for  a  two-dimensional  flat  plate. 

The  theory  of  Bollay  is  limited  to  rectangular  wings  of  small 
aspect  ratio  in  incompressible  flow  and  is  not  suited  to  the  calcu¬ 
lation  of  quantities  other  than  normal  force.  However,  the  agree¬ 
ment  between  the  calculated  and  measured  normal  forces  is  excellent 
at  very  low  aspect  ratios  up  to  angles  of  attack  of  about  40°. 
Furthermore ,  the  predicted  normal  forces  for  high  aspect  ratios 
(A  ~  6)  agree  with  experiment  within  about  20  percent  for  angles 
of  attack  below  the  stall.  The  stalling  phenomenon  for  high  aspect 
ratios  is  clearly  contrary  to  the  assumed  mathematical  model. 

A  few  years  after  the  appearance  of  Bollay 's  paper,  a  theo¬ 
retical  analysis  was  published  by  F.  Weinig  (Ref.  3)  whose  aim  was 
to  extend  the  Prandtl  wing  theory  down  to  lower  aspect  ratios. 

Weinig  observed  that  for  small-aspect-ratio  wings  the  apparent 
mass  of  air  which  is  deflected  downward  by  the  wing  is  greater  than 
that  which  results  from  lifting  line  theory  with  elliptic  lift  dis¬ 
tribution.  According  to  lifting  line  theory,  this  mass  would  corre¬ 
spond  to  the  mass  of  a  cylinder  of  air  having  a  diameter  equal  to 
the  span.  Weinig 's  innovation  was  twofold:  first,  he  considered 
the  actual  projected  area  of  the  wing  (S  sin  a)  so  that  the 
mass  of  air  deflected  downward  per  unit  time  is  given  by 

m  ■  pV  (irs^  +  S  sin  a)  (5) 

*  00 

where  s  is  the  semispan  and  S  the  planform  area;  and  second, 
he  introduced  a  novel  method,  based  on  cascade  theory,  for  calcu¬ 
lating  the  induced  downwash  associated  with  this  mass . 

In  this  manner,  Weinig  developed  a  theory  for  the  lift  and 
drag  of  small -aspect-ratio  wings  which  is  considerably  simpler 
than  Bollay 's  theory  from  a  mathematical  standpoint  and  which  gives 
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identical  results  (Eqs.  (3)  and  (4))  for  the  limiting  cases  of 
infinite  and  zero  aspect  ratio.  The  agreement  with  experiment 
for  nonvanishing  aspect  ratios,  however,  is  not  quite  as  good  as 
for  Bo  Hay's  theory. 

2.1.2  Intermediate  aspect  ratio 

More  recently,  K.  Gersten  (Ref.  4)  presented  an  analytical 
method  for  extending  linear  lifting  surface  theory  down  to  lower 
aspect  ratios  by  incorporating  the  low  aspect  ratio  results  of 
Bollay.  That  is,  Gersten  allows  the  trailing  vortices  of  each 
elementary  horseshoe  vortex  to  be  shed  above  the  plane  of  the  wing. 
Thus,  the  actual  wing  is  considered  as  made  up  of  a  number  of 
Bollay-type  wings  placed  side  by  side  and  it  is  assumed  that  the 
shedding  angle  of  all  the  horseshoe  vortices  remains  constant  at 
the  limiting  value  of  a/2  found  by  Bollay  for  vanishing  aspect 
ratios  and  angles  of  attack.  Gersten* s  model  therefore  has 
vortices  shed  over  the  entire  wing  surface.^  In  order  to  make  the 
problem  mathematically  tractable,  Gersten  then  divided  the  wing 
into  a  finite  number  of  spanwise  strips,  concentrated  the  lifting 
line  of  each  strip  at  the  quarter  chord,  and  proceeded  to  satisfy 
the  boundary  condition  of  no  flow  through  the  wing  at  the  corre¬ 
sponding  three-quarter  chord  of  each  strip.  Thus,  Gersten 's 
mathematical  model  consists  of  a  discrete  number  of  lifting  lines, 
each  having  an  unspecified  spanwise  loading  and  shedding  a  flat, 
continuous  vortex  sheet  at  half  the  angle  of  attack,  as  shown 
below; 


^This  type  of  vortex  shedding  is  incompatible  with  the  physical 
requirement  that  the  vortices  follow  the  local  streamlines,  which 
must  lie  in  the  surface  of  the  wing.  (See  discussion  on  page  16.) 
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Y_^(n) 


Mathematically,  then,  Gersten  represents  the  spanwise  circu¬ 
lation  distribution  of  each  lifting  line  as  a  sine  series  and,  by 
applying  the  boundary  condition  at  discrete  points  on  each  3/4- 
chord  line,  is  led  to  a  set  of  simultaneous  algebraic  equations 
for  the  unknown  coefficients.  In  arriving  at  these  equations, 
Gersten  further  assumed  that  the  angle  of  attack  is  small  and 
neglected  higher  order  terms  in  angle  of  attack.  Then,  the  normal 
force  and  the  pitching  moment  can  be  represented  as 

N  ■■  N  +  N 

linear  theory  sep 

-  k  a  +  k  (6) 

1  s 


and 


M  ■  k  a  +  k 

3  4 


(7) 
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so  that  the  solution  of  Reference  4  is  actually  made  up  of  the 
linear  theory  of  Scholz  (Ref.  5)  plus  a  quadratic  correction  term 
which  is  determined  from  the  analysis  and  stems  from  the  fact  that 
the  trailing  vortices  do  not  lie  in  the  plane  of  the  wing. 

It  is  clear  that  the  above  model  does  not  reduce  to  Bollay's 
model  for  very  low  aspect  ratios.  In  fact  at  aspect  ratio  zero, 
Gersten's  model  yields  a  coefficient  of  3.2  for  the  quadratic  term 
in  the  normal  force,  as  compared  with  Bollay's  value  of  2.0. 
Gersten's  mathematical  procedure  is  by  nature  restricted  to  small 
angles  of  attack,  but  appears  to  give  an  improvement  over  Bollay's 
theory  for  the  aspect  ratio  range  from  2  to  4  for  angles  of  attack 
up  to  about  15°.  In  addition,  Gersten's  analysis  is  suited  to  the 
calculation  of  load  distribution  and  pitching  moment  in  addition 
to  lift  and  drag. 

2.1.3  High  aspect  ratio 

For  aspect  ratios  of  the  order  of  5  or  6  and  above,  with 
rounded  wing  tips,  the  flow  over  the  wing  remains  attached  up  to 
the  stall.  Below  the  stall,  the  classical  lifting  line  theory 
applies  and  the  lift  curve  is  essentially  linear.  The  stall 
itself  becomes  more  abrupt  as  the  aspect  ratio  increases  and  does 

notVbehave  according  to  the  steady-flow  models  proposed  by  Bollay 

/■ 

and  Gersten.  Instead,  additional  vortices  which  are  aligned  normal 
to  the  stream  direction  begin  shedding  from  the  leading  and  trail¬ 
ing  edges,  and  the  flow  becomes  unsteady.  The  problem  then  becomes 
the  bluff  body  problem,  which  remains  one  of  the  classical  unsolved 
problems  of  aerodynamics.  There  is,  therefore,  as  yet  no  satis¬ 
factory  theory  for  treating  high-aspect-ratio  wings  above  the  stall 
angle,  which  is  well  defined  experimentally.  However,  there  has 
been  some  interesting  analytical  work  related  to  this  problem 
(Refs.  6  and  7)  which  employs  the  classical  integral-differential 
equation  of  Prandtl  lifting-line  theory.  The  essential  contri¬ 
bution  of  this  work  is  a  demonstration  that  if  the  airfoil  section 
is  assumed  to  have  a  lift  curve  which  shows  discontinuities  or 
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negative  slopes  (which  are  typical  of  the  stall) ,  then  the  Prandtl 
equation  yields  asymmetric  solutions  for  the  span  loading  and 
therefore  predicts  the  possibility  of  large  rolling  moments  asso¬ 
ciated  with  the  stall.  The  thesis  work  of  Schairer  (Ref.  7)  treats 
a  specific  case  of  a  discontinuous  lift  curve,  as  shown  in  the 
following  sketch; 


2.2  Triangular  Wings 

The  theory  of  Brown  and  Michael  (Ref.  8)  contains  the  essential 
elements  of  a  number  of  recent  papers  (e.g.,  Refs.  9,  10,  and  11) 
which  collectively  represent  the  present  status  of  the  mathematical 
analysis  of  highly  swept  wings  with  leading-edge  separation.  The 
analysis  is  confined  to  triangular  wings  of  low  aspect  ratio  and 
makes  use  of  the  s lender-body  approximation.  With  this  approxi¬ 
mation,  the  total  lift  is  given  by  (see  Ref.  12) 

L  -  -pV„  j  <p  dy  (8) 

c 

where  the  contour  c  encloses  the  trailing-edge  cross  section  of 
the  wing  and  the  vortex  sheets.  Thus,  the  lift  is  linear  in  the 
potential  0  and  is  therefore  composed  of  two  parts;  one  being 
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the  usual  slender-body  lift  (which  is  linear  in  a) ,  and  one  being 
the  additional  lift  due  to  the  potential  associated  with  the  sepa¬ 
rated  vortices. 

The  model  proposed  by  Brown  and  Michael  is  illustrated  in  the 
following  sketch: 


It  is  assumed  that  all  of  the  trailing  vorticity  above  the  wing  is 
concentrated  in  two  fully  rolled-up  vortices  whose  strengths  vary 
linearly  in  the  chordwise  direction.  The  contribution  of  the  feed¬ 
ing  vortex  sheets  (between  the  wing  leading  edges  and  the  rolled-up 
vortices)  to  the  co’^plex  potential  in  the  cross-flow  plane  is 
neglected.  However,  since  the  feeding  vortices  shed  from  the  lead¬ 
ing  edge  lie  almost  normal  to  the  free  stream  (rather  than  along 
the  streamlines)  they  must  sustain  a  force.  The  condition  that  the 
shed  vortex  system  be  force -free  everywhere  is  therefore  approxi¬ 
mated  by  requiring  that  the  net  force  on  the  entire  vortex  system 
be  zero.  In  this  manner,  with  the  Kutta  condition  imposed  at  the 
leading  edges.  Brown  and  Michael  proceed  to  solve  for  the  vortex 
strengths  r  and  the  vortex  positions,  and  thereby  calculate  the 
total  lift  and  drag,  as  well  as  the  spanwise  load  distribution^. 


similar  analysis  has  been  applied  by  Cheng  (Ref.  13)  to  the 
corresponding  problem  of  a  rectangular  wing  with  side-edge 
separation  in  supersonic  flow.  Here  the  separation  effects  are 
confined  to  the  tip  region  within  the  Mach  cone. 
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The  lift  curve  calculated  by  the  above  method  agrees  well 
with  experiment  only  for  extremely  low-aspect  ratios  (small  apex 
angles)  and  the  predicted  high  peaks  in  the  pressure  distribution 
are  not  realized  because  the  streamwise  vorticity  is  actually 
distributed  through  the  feeding  sheet  rather  than  being  concen¬ 
trated  in  two  discrete  vortices. 

A  somewhat  more  general  treatment  of  the  slender  triangular 
wing  with  leading -edge  separation  has  been  developed  by  Mangier 
and  Smith  (Ref.  14)  who  consider  a  more  realistic  shape  of  the 
vortex  sheets  and  calculate  the  lateral  distribution  of  streamwise 
circulation  within  the  sheet,  accounting  approximately  for  the 
requirement,  of  no  force  on  the  sheet  and  no  flow  through  it.  (See 
sketch . ) 


However,  the  severe  limitation  on  aspect  ratio  remains. 

2 . 3  Sweptback  Wings 

In  Reference  15,  Pappas  and  Kunen  treated  the  high-aspect- 
ratio  swept  wing  by  what  is  essentially  simple  sweep  theory  with 
a  single  concentrated  vortex  representing  the  leading-edge  separa¬ 
tion.  Thus,  the  flow  in  planes  normal  to  the  leading  edge  is 
treated  as  shown  in  the  sketch  below. 
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and  the  lift  and  drag  of  the  airfoil  section  are  calculated  by  the 
two-dimensional  Blasius  theorem  involving  the  complex  potential 
in  the  plane  A-A.  This  treatment  is  similar  to  the  s lender-body 
analysis  of  Brown  and  Michael  (Ref.  8)  for  the  slender  delta  wing 
where  the  two-dimensional  flow  in  the  cross-flow  plane  is  calcu¬ 
lated.  Hence,  there  are  similar  dravrbacks  to  the  analysis  of  Pappas 
and  Kunen  (Ref.  15)  who  calculate  the  separation  vortex  strength 
r  from  the  experimentally  determined  vortex  position  over  the  swept 
wing.  That  is,  the  separated  vortex  strength  r  must  actually 
vary  spanwise  over  the  sweptback  wing.  Therefore,  there  must  be 
a  feeding  vortex  sheet  connecting  the  separated  vortex  with  the 
wing  leading-edge,  and  the  force  acting  on  this  sheet  plus  the 
separated  vortex  must  be  set  to  zero.  This  was  not  done  in  Ref¬ 
erence  15.  However,  good  agreement  with  experiment  was  demonstrated 
for  the  calculated  span  load  distribution  over  the  inner  half  of 
the  span  for  the  case  of  a  45°  swept  wing  at  10°  angle  of  attack. 

The  theory  of  Reference  15  is  clearly  restricted  to  high 
aspect  ratios  and  does  not  apply  over  the  outboard  sections  of 
the  wing,  where  the  flow  pattern  is  quite  different  from  the 
assumed  two-dimensional  model.  Pappas  and  Kunen  suggest  using 
experimental  stalled-airfoil  data  in  the  tip  regions. 
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Just  in  the  past  year,  an  attempt  has  been  made  to  treat 
sweptback  wings  of  arbitrary  aspect  ratio  with  flow  separation. 

This  analysis  (Ref.  16)  by  Gersten  represents  a  natural  extension 
of  his  rectangular  wing  analysis  (Ref.  4)  which  was  discussed  in 
Section  2.1.2.  In  Reference  16,  however,  the  Bollay  model  is  com¬ 
bined  with  the  lifting  surface  model  of  Truckenbrodt  (Ref.  17)  to 
produce  a  wing  element  with  side-edge  separation.  That  is,  the 
trailing  portion  of  the  horseshoe  vortices  making  up  each  chord- 
wise  strip  in  the  lifting  surface  theory  are  permitted  to  drift 
above  the  wing  surface  at  half  the  angle  of  attack.  Thus,  Gersten 
again  permits  vortex  shedding  over  the  entire  wing  surface.  In 
Reference  16,  as  in  Reference  4,  the  normal  force  is  expressed  as 
a  linear  term  plus  a  quadratic  term,  the  former  now  being  given 
by  the  lifting  surface  theory  of  Truckenbrodt.  The  chordwise 
loading  is  assumed  to  be  given  by  the  first  two  terms  of  the 
Birnbaum  series,  and  the  boundary  condition  of  no  flow  through 
the  wing  is  satisfied  at  both  the  quarter-chord  and  the  wing 
trailing  edge. 

In  this  manner,  Gersten  has  built  up  a  model  of  arbitrary 
planform  and  aspect  ratio  with  a  solid  wake  of  streamwise  vorticity 
forming  an  angle  a/2  above  the  plane  of  the  wing  (see  sketch) . 


The  agreement  between  theory  and  experiment  for  Gersten 's  model 
appears  to  be  satisfactory,  although  better  agreement  for  swept 
and  delta  wings  of  aspect  ratios  below  unity  is  obtained  with 
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Newtonian  impact  theory*,  and  still  better  agreement  for  very  low 
aspect  ratio  rectangular  wings  is  obtained  with  the  theory  of 
Bollay  (Ref.  1) . 

2.4  Summary  of  State  of  the  Art 

From  the  foregoing  review  of  the  analytical  papers  pertinent 
to  the  problem  of  wings  with  separated  flow,  it  is  apparent  that 
the  theory  is  lacking  in  several  respects .  With  the  exception  of 
Reference  15  (which  does  not  give  total  forces)  and  Reference  16 
(which  violates  a  fundamental  boundary  condition) ,  the  analytical 
work  on  sweptback  wings  with  leading-edge  separation  has  been 
confined  to  delta  wings  of  extremely  low  aspect  ratio  and  shows 
little  improvement  over  the  original  s lender-wing  analysis  of 
R.  T.  Jones  (Ref.  18)  insofar  as  the  prediction  of  total  aerody¬ 
namic  forces  is  concerned,  except  for  aspect  ratios  less  than  one. 
These  theories  do,  however,  give  a  more  realistic  prediction  of 
the  span  loading  on  a  very  slender  delta  wing  at  moderate  angles 
of  attack,  indicating  the  peaks  under  the  separated  vortices,  as 
opposed  to  the  elliptic  span  load  distribution  of  the  attached- 
flow  solution  of  Jones. 

For  the  rectangular  wing,  the  theory  of  Bollay  is  similarly 
restricted  to  rather  low  aspect  ratios  but  is  valid  up  to  high 
angles  of  attack  and  shows  excellent  agreement  with  the  experi¬ 
mentally  measured  normal  force  on  a  rectangular  wing  of  aspect 
ratio  A  “  1/30  up  to  an  angle  of  attack  of  about  40°.  If 
Bollay' s  theory  is  to  be  properly  extended  to  higher  aspect  ratios, 
however,  one  must  clearly  abandon  both  the  assumption  of  uniform 
span  loading  and  that  of  the  two-dimensional  form  of  the  chord- 
wise  loading.  These  steps  were  both  taken  by  Gersten  in  Refer¬ 
ences  4  and  16,  but  the  resulting  mathematical  model  violates 
certain  physical  requirements  which  are  fundamental  to  the  validity 


*Vidya  is  publishing  a  separate  report  on  this  subject,  Vidya 
Report  No.  92  by  Sacks  and  Burnell,  dated  March  30,  1963. 
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of  the  results.  The  agreenent  of  Gersten's  theory  with  experiment, 
which  extends  only  to  moderate  angles  of  attack,  must  therefore 
be  considered  fortuitous. 

The  objections  to  Gersten's  model  can  most  easily  be  under¬ 
stood  if  we  consider  the  physical  problem  in  somewhat  greater 
detail.  Perhaps  the  most  basic  observation  is  that  there  can  be 
no  flow  of  fluid  through  a  solid  surface.  Thus,  a  boundary  con¬ 
dition  on  the  mathematical  problem  is  that  the  normal  component 
of  velocity  at  the  wing  surface  must  be  zero®.  In  other  words, 
the  wing  is  a  stream  surface  and  the  streamlines  must  lie  in  the 
surface  of  the  wing  except  at  points  of  flow  separation.  Further, 
according  to  Helmholtz's  law  of  vortex  motion,  a  free  vortex  must 
follow  the  streamlines.  Therefore,  vortices  can  be  shed  only 
along  lines  of  flow  separation.  For  the  sharp-edged  flat  plate, 
the  lines  of  flow  separation  coincide  with  the  edges  of  the  plate, 
and  vortex  shedding  can  therefore  occur  only  along  the  sharp 
edges  1  It  can  be  seen  that  the  mathematical  model  postulated  by 
Gersten  implies  an  arbitrary  number  of  spanwise  lines  cf  flow 
separation  over  the  upper  surface  of  the  wing  and  is  evidently 
unacceptable.  Bollay's  model,  on  the  other  hand,  permits  flow 
separation  only  along  the  side  edges.  It  should  also  be  pointed 
out  that  in  both  of  these  models  it  has  been  assumed  that  all  of 
the  vorticity  is  shed  out  of  the  plane  of  the  wing.  This  means 
that  there  are  no  chordwise  vortices  in  the  plane  of  the  wing  and 
hence  there  is  no  vortex  sheet  lying  in  the  wing  and  shed  at  the 
trailing  edge!  Because  of  this  assumption,  it  was  not  necessary 
in  the  analyses  of  References  1,  4,  ai.J  16,  to  in^ose  the  Kutta 
condition  along  the  side  edges,  since  the  velocities  there  were 
necessarily  finite.  However,  if  chordwise  vortices  in  the  plane 


®This  boundary  condition  is  ordinarily  satisfied  only  at  a  finite 
nunber  of  points  on  the  wing  surface  (as  in  Ref.  4)  or  else  in 
the  mean  (as  in  Ref.  1). 
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of  the  wing  are  permitted,  as  they  should  be,  then  the  Kutta  condi¬ 
tion  must  be  imposed  along  the  side  edges  to  insure  that  the 
velocities  there  remain  finite.  This  is  equivalent  to  the  physical 
requirement  that  the  flow  must  separate  tangentially  at  all  sharp 
edges . 

In  the  following  sections,  a  new  theory  will  be  developed  for 
rectangular  wings  with  side-edge  separation  which  will  in  fact 
satisfy  these  fundamental  requirements  and  therefore  will  furnish 
a  sound  basis  upon  which  we  shall  be  able  to  develop  a  treatment 
for  other  planforms  as  well.  The  flow  will  be  considered  to  be 
three-dimensional,  so  that  no  restriction  on  aspect  ratio  is  implied. 

3.  APPROACH  TO  THE  PROBLEM 

The  problem  to  be  treated  in  the  present  analysis  is  that  of 
a  sweptback  wing  of  arbitrary  aspect  ratio  in  incompressible  flow 
with  separation  along  the  entire  leading  edge.  Such  flow  separa¬ 
tion  is  steady  up  to  large  angles  of  attack,  depending  upon  the 
sweep  angle,  but  is  three-dimensional  in  character,  becoming  conical 
only  for  sweep  angles  near  90°  (i.e.,  for  vanishing  aspect  ratios). 
The  use  of  strip  theory  in  either  the  chordwise  or  spanwise  direc¬ 
tion  is  therefore  not  appropriate.  However,  if  the  trailing  edge 
is  straight,  then  we  can  represent  a  sweptback  wing  of  arbitrary 
aspect  ratio  and  rather  general  planform  by  a  number  of  rectangular 
wing  elements  (see  sketch) ,  provided  that  each  element  satisfies 
the  appropriate  boundary  and  side-edge  conditions. 
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Since  separation  from  the  equivalent  wing  (composed  of  rec¬ 
tangular  elements)  can  occur  only  along  its  outline  (i.e.,  along 
its  exposed  edges),  leading-edge  separation  can  be  permitted  only 
over  the  exposed  portion  of  the  span  of  each  rectangular  wing. 
Furthermore,  leading -edge  separation  from  a  rectangular  wing  occurs 
at  higher  angles  of  attack  and  is  essentially  an  unsteady  phenomenon 
in  which  the  separated  vortices  are  normal  to  the  free  stream  and 
are  shed  periodically.  Therefore,  for  the  purpose  of  treating  the 
present  problem,  flow  separation  will  be  permitted  only  along  the 
side  edges  of  each  rectangular  wing  element.  Thus,  as  the  number 
of  wing  elements  is  increased,  the  mathematical  approximation  will 
approach  the  physical  problem  of  steady  flow  separation  along  a 
sweptback  leading  edge  (see  sketch) . 


In  the  present  analysis,  each  rectangular  wing  element  will 
be  treated  in  the  manner  of  classical  lifting  line  theory,  except 
that  side-edge  flow  separation  will  be  included  and  the  Kutta 
condition  will  be  satisfied  along  each  side  edge.  This  latter 
point  is  essential  to  the  success  of  the  basic  approach,  since  it 
automatically  removes  the  chordwise  singularity  in  the  downwash 
along  each  side  edge,  which  arises  in  Prandtl  lifting  line  theory. 
Only  by  the  removal  of  these  singularities  can  one  safely  super¬ 
pose  rectangular  elements  in  the  manner  sketched  above. 
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It  should  be  noted  that,  since  the  rectangular  wing  is  repre¬ 
sented  by  a  single  lifting  line  (plus  a  separated  vortex  system) , 
one  cannot  expect  an  accurate  prediction  of  the  chordwise  load 
distribution  or  of  the  center  of  pressure.  However,  these  details 
can  be  obtained  if  desired  by  building  up  a  rectangular  wing  with 
a  nuitber  of  higher  aspect  ratio  rectangular  elements  in  the  same 
manner  as  for  the  swept  wing  (see  sketch  below) . 


In  the  following  section,  the  theoretical  analysis  will 
first  be  developed  for  the  basic  rectangular  wing  with  side -edge 
flow  separation.  This  analysis  will  then  be  applied  to  swept- 
back  wings  with  straight  trailing  edges  by  superposition,  using 
the  rectangular  wing  as  an  elementary  building  block  in  the 
manner  described  above.  In  addition,  the  extension  of  the  analy¬ 
sis  to  incorporate  the  influence  of  the  ground  plane  will  be 
outlined  in  an  appendix. 
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4.  THEORY  OF  RECTANGULAR  WINGS  WITH  SIDE-EDGE  SEPARATION 
4.1  Statement  and  Discussion  of  the  Problem 

The  problem  to  be  treated  in  this  section  is  that  of  a  rec¬ 
tangular  wing  element  of  arbitrary  aspect  ratio  in  incompressible 
flow  with  side-edge  separation.  Separation  from  the  leading  edge 
will  not  be  considered,  since  it  is  intended  to  apply  the  theory 
later  to  sweptback  wings,  as  explained  in  the  previous  paragraphs. 
Thus,  the  rectangular  wing  will  be  represented  mathematically  by 
(1)  a  lifting  line  of  unspecifjd  spanwise  circulation  distribu¬ 
tion  7  (ti)  ,  (2)  its  associated  trailing  vortex  sheet  of  strength 

-(dy/dr|)  dii  lying  in  the  plane  of  the  wing,  and  (3)  a  continuous 
system  of  horseshoe  vortices  of  strength  (df/d^)  d|  representing 
the  separated  vortex  system.  This  representation  is  shown  sche¬ 
matically  in  Figure  1. 

It  is  assumed  that  separation  occurs  all  along  the  side 
edges  and  that  all  of  the  vortices  shed  from  the  side  edges  lie 
in  the  two  planes  containing  the  free-stream  direction  and  the 
side  edges,  in  the  manner  of  Bollay.  It  is  evident  that  the 
boundary  conditions  of  no  flow  through  the  wing  can  be  satisfied 
along  a  selected  control  line  by  a  suitable  choice  of  7(11)  for 
any  given  dF/d^.  Therefore,  the  Kutta  condition  of  finite  veloc¬ 
ity  at  the  side  edges  will  have  to  be  imposed  to  render  dF/d^ 
(which  will  be  assumed  constant  in  the  present  analysis)  unique. 
Here  again,  the  justification  for  this  assumption  stems  from  the 
intended  application  to  swept  wings  wherein  dF/d^  may  vary  from 
one  rectangular  wing  element  to  the  next.  Hence,  the  assumption 
that  dF/d^  is  constant  on  a  given  rectangular  wing  element  is 
of  no  consequence  in  the  sw$pt  wing  problem . 

Again,  it  is  pointed  out  that  the  representation  of  the  rec¬ 
tangular  wing  by  a  single  lifting  line  (plus  a  separated  vortex 
system  of  constant  strength  per  unit  length)  may  later  be  improved 
by  the  use  of  a  number  of  such  rectangular  wing  elements  of  higher 
aspect  ratio,  placed  one  behind  the  other. 
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Note  that  the  above  mathematical  model  satisfies  the  physical 
conditions  discussed  in  Section  2.4  and  differs  from  the  models  of 
Bollay  (Ref.  1)  and/or  Gersten  (Ref.  4)  in  the  following  in^rtant 
respects : 

(1)  The  span  load  distribution  is  not  specified. 

(2)  The  appropriate  streamwise  vorticity  in  the  plane  of 

the  wing  is  included.  Thus  there  is  a  vortex  sheet 
shed  from  the  trailing  edge. 

(3)  Vortex  shedding  occurs  only  at  the  side  and  trailing 

edges  of  the  wing . 

4.2  Analysis 

4.2.1  Mathematical  formulation  of  boundary-valutf  problem 

The  mathematical  equations  to  be  solved  for  the  above  problem 
are  actually  those  which  express  the  boundary  condition  of  no  flow 
through  the  wing  at  each  of  a  number  of  selected  control  points 
on  the  wing  surface.  Thus,  if  we  express  the  upwash  due  to  the 
lifting  line  and  its  trailing  sheet  in  the  plane  of  the  wing  as 
w*(x,y)  and  the  upwash  due  to  the  horseshoe  vortices  of  the  sepa¬ 
rated  system  as  w(x,y) ,  then  the  boundary  condition  at  a  selected 
point  wing  becomes 

w*(Xp,yp)  +  w(Xp,yp)  -  -V^  sin  a  (9) 

Ordinarily,  this  tangency  condition  of  no  flow  through  the 
wing  is  satisfied  at  a  number  of  discrete  points,  and  the  number 
of  simultaneous  equations  to  be  solved  is  then  equal  to  the  number 
of  points  selected.  However,  in  the  present  analysis,  we  shall 
satisfy  Equation  (9)  for  each  Fourier  component.  Thus,  the  number 
of  simultaneous  equations  to  be  solved  is  equal  to  the  number  of 
terms  tedeen  in  the  Fourier  analysis,  and  the  boundary  condition  is 
satisfied,  approximately,  for  all  spanwise  points  along  the  selected 
control  line®. 


®This  method  avoids  the  question  of  where  control  points  should 
be  located.  Also  its  accuracy  can  be  assessed  by  examining  the 
size  of  subsequent  Fourier  terms. 
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The  following  analysis  will  therefore  consist  of  setting  up 
mathematical  expressions  for  w*(x,y)  and  ^(x,y),  performing  a 
Fourier  analysis  in  a  spanwise  variable  related  to  y,  and  sub¬ 
stituting  these  e}q>ressions  into  Equation  (9),  which  will  then  be 
satisfied  for  each  Fourier  con^nent.  Finally,  the  resulting  set 
of  simultaneous  equations  will  be  solved  for  the  unknown  Fourier 
coefficients  of  the  expansion  representing  the  spanwise  circula¬ 
tion  distribution  7(11)*  Once  this  loading  is  obtained,  the 
strength  of  the  separated  vortex  system  will  be  determined  uniquely 
by  imposing  the  Kutta  condition  along  the  side  edges  y  ~  ±  s, 
and  the  resulting  aerodynamic  forces  and  moments  will  be  calculated. 

4.2.2  Velocity  field  of  separated  horseshoe  vortex  system 

Each  horseshoe  vortex  making  up  the  separated  vortex  system 
is  composed  of  three  parts,  a  bound  vortex  in  the  wing  and  a  left 
and  a  right  trailing  vortex  lying  at  an  angle  6  above  the  wing. 

It  is  therefore  convenient  to  write  the  corresponding  upwash  in  the 
plane  of  the  wing  w(x,y)  in  three  parts  as 

W(x,y)  -  Wg(x,y)  +  g^(x,y)  +  Wg(x,y)  (10) 

The  expression  for  each  of  these  con^nents  at  an  arbitrary  point 
in  space  has  been  derived  in  Appendix  A  for  a  single  horseshoe 
vortex  of  strength  F.  Thus,  replacing  F  in  the  expressions 
of  Appendix  A  by  (dF/d|)  d^,  setting  z  to  zero,  and  integrating 
over  the  chord  c,  we  have  for  the  continuous  vortex  distribution 
representing  the  side-edge  separation  (see  Fig.  1) 
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In  the  present  analysis,  we  shall  assume  that  dr/d^  is  con¬ 
stant  along  the  chord,  which  means  that  vortex  shedding  is  uniform 
along  the  side  edges.  (In  the  case  of  the  swept  wing  to  be  treated 
in  Section  5,  the  values  of  6  and  dr/d^  will  be  permitted  to 
vary  from  one  rectangular  wing  element  to  the  next.)  Thus,  noting 
that  the  integral  for  Wg(x,y)  is  singular  for  0  <  x  <  c  and 
taking  the  Cauchy  principal  value  for  that  integral,  we  find  upon 
integration 


(12a) 


These  three  expressions,  then,  represent  the  upwash  at  any  point 
(x,y)  in  the  plane  of  the  wing  due  to  the  separated  vortex  system 
shown  in  Figure  1. 

4.2.3  Velocity  field  of  vortex  system  in  wing  surface 

It  will  be  convenient,  as  in  classical  lifting  line  theory 
(e.g..  Ref.  19)  to  represent  the  unknown  circulation  distribution 
^(t))  along  the  lifting  line  as  a  Fourier  sine  series  of  the  form 

(13) 


(14) 


'y(-V')  ■  -yif) 


(15) 
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Thus,  the  wing  lies  in  the  interval  -s  <  tj  ^  s  which  corresponds 
to  0  ^  f  ^  TT .  It  can  be  seen  that  the  series  of  Equation  (13) 
satisfies  the  requirement  that  y  be  zero  at  the  wing  tips,  regard¬ 
less  of  how  few  terms  are  used  in  the  series,  and  that  the  selected 
series  is  even  in  t],  which  insures  lateral  symmetry. 

If  we  consider  now  that  the  circulation  distribution  yCt)) 
can  be  thought  of  as  a  large  number  of  steps  as  shown  in  the  sketch. 


-S  +S  T) 

0  ir  f 


then  it  is  clear  that  a  vortex  of  strength  -(ciY/dr|)  dt)  must  be 
shed  at  each  step  and  we  can  therefore  construct  the  entire  vortex 
system  in  the  plane  of  the  wing  from  elementary  horseshoe  vortices 
of  semispan  s^  and  strength  -(dry/dr))  dr]  as  shown  below: 


z 
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Thus,  if  we  represent  the  upwash  (in  the  plane  of  the  wing)  due 
to  an  elementary  horseshoe  vortex  of  unit  strength  as  w^,  we 
have  for  the  total  upwash  at  a  point  (x,y)  on  the  wing  due  to 
the  entire  vortex  system  in  the  plane  of  the  wing 

s 

w*(x,y)  ”  -  J  w^(x,y,^i,s^)  ^  di^  (16) 

o 

Upon  introducing  Equation  (13)  for  7(11)  and  interchanging  the 
order  of  integration  and  summation,  we  have 
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Now  we  can  employ  the  expressions  of  Appendix  A  for  the  upwash 
w^  due  to  each  elementary  horseshoe  vortex  by  setting  z  ■  0 
and  letting 

r  -  1 
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Thus  we  find  from  Equations  (A. 15)  that 
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The  total  upwash  w*  due  to  Y(n)  is  therefore  given  by 
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Now  it  can  be  shown  by  a  simple  change  of  variable  and  algebraic 
manipulation  that  for  odd  values  of  n 
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80  that  Equation  (20)  can  be  collapsed  into  two  integrals  over  the 
larger  interval  of  0  to  tt.  Thus,  conibining  the  four  terms 
involving  the  square  root,  and  making  use  of  Equations  (21),  we  find 
that  Equation  (20)  can  be  written  in  the  following  form: 
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Both  of  the  integrals  in  Equation  (22)  are  improper  (i.e., 
singular)  integrals.  However,  i'f  we  replace  the  control-point 
variable  y  by  (-s  cos  V'q)  »  then  we  recognize  the  first  integral 
as  that  of  classical  lifting-line  theory.  For  points  within  the 
span,  the  value  of  this  integral  is  (Ref.  19) 
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The  singularity  of  the  second  integral  in  Equation  (22)  can  be 
removed  by  sinply  adding  and  subtracting  a  term.  Thus 


/ 


J  V(x-|j^)^-»-(y+s  cos  f)' 


(x-|^) (y+s  cos  f) 


cos  nf  df 


J  I  (x-^^^)(y+8  cos  V') 


cos  n^  df/f 


TT 

‘J 


Si 

y*8  cos  f 


(24) 


-29- 


where  the  lower  sign  applies  if  <  0  (i.e.,  for  points 

ahead  of  the  lifting  line)  and  we  recognize  that  the  last  Integral 
is  identical  with  that  of  Equation  (23).  The  remaining  integral 
is  non-singular  and  is  furthermore  recognizable  as  the  integral 
giving  the  Fourier  coefficients  of  a  certain  cosine  series  in 
That  is,  if  we  let 
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where  the  lower  sign  again  applies  if  x  -  <  0,  then 
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and  the  resulting  relation  is 


r  cos  V')®  ^  /x-^i  \ 

J  “Tx-'epTy+s  cos  f) -  nf  df  >  ^oj  ^ 


TT  sin  nf^ 
s  sin  f 


(27) 


Substitution  of  Equations  (23)  and  (27)  into  Equation  (22)  therefore 
yields  the  final  expression  for  w*  (for  points  within  the  span)  in 
the  form 
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The  first  term  of  Equation  (28a)  is  that  appertaining  to  a  value 
of  (x  -  4j^)/s  which  is  sufficiently  far  behind  the  lifting  line 
that  the  upwash  no  longer  depends  on  (x  -  ^^)/s.  This  is  evidently 
twice  the  upwash  at  the  lifting  line  itself. 

It  is  instructive  to  look  at  the  term  n  -  1  of  the  foregoing 
equations.  The  circulation  distribution  from  Equations  (13)  and 
(14)  is  then 
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so  that  the  corresponding  span  loading  is  elliptical.  Thus  letting 


f 


n 


^i 


27r  sin 

_ la 

sin 

o 


+ 


]L 

2 


(30) 


we  find  that  the  upwash  behind  the  lifting  line  is 
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Table  I  shows  the  variation  of  downwash  (-w*/V^)  with  longitudinal 
distance  behind  the  lifing  line  in  semispans  and  with  lateral 
distance  from  the  wake  centerline  in  semispans.  It  is  interesting 
to  note  that  one  semispan  behind  the  lifting  line,  the  downwash 
is  still  about  10  percent  higher  than  its  value  infinitely  far 
downstream  (where  f^  >  2ir)  and  that  it  is  quite  uniform  across 
the  span. 
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4.2.4  Solution  of  the  boundary-value  problem 

The  set  of  equations  representing  the  present  boundary  value 
problem  singly  express  the  boundary  condition  that  there  can  be 
no  flow  through  the  wing  surface.  Thus,  having  derived  the  mathe¬ 
matical  expression  for  the  normal  velocity  w  due  to  the  separation 
vortices  and  w*  due  to  the  lifting  line  and  its  trailing  sheet, 
we  can  express  this  condition  as 

^  ^  ^*p*yp^  -  -  sin  a  (32) 

That  is,  the  velocity  component  normal  to  the  wing  surface  due  to 
all  vortices  must  exactly  cancel  the  normal  component  of  the  free- 
stream  velocity.  The  point  point  on  the  wing  at 

which  this  condition  is  to  be  satisfied,  and  the  number  of  equations 
to  be  solved  for  the  unknown  distribution  'y(Ti)  is  ordinarily 
equal  to  the  selected  number  of  control  points.  However,  in  the 
present  analysis,  we  shall  eliminate  the  dependence  on  y^  by 
Fourier  analysis  and  thus  satisfy  Equation  (32)  for  all  y^  within 
the  span  at  a  specified  control  line  x  ■  x^.  The  number  of  equa¬ 
tions  will  therefore  be  equal  to  the  number  of  Fourier  coicqponents 
employed. 

In  order  to  expand  w*  in  a  Fourier  series  of  it  is 

first  convenient  to  note  from  the  trigonometric  identity  (Ref.  20) 

sin  mx  cos  nx  ■  j  sin(n  +  m)x  -  ^  sin(n  -  m)x 
that  sin  nf^/sin  can  be  written  in  the  form 


sin  nV'jj 
sin 

’O 


n-i 
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>2,4 


COS  jfjj 
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“IT  1  V'o  1  IT 


1  +  2 


(33) 
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Therefore,  we  shall  expand  w*  and  w  in  the  form  of  cosine 
series  in  Furthermore,  only  even  multiples  of  can 

appear  because  of  the  required  symmetry  about  the  midspan 

-  (f/2)  .  Thus,  we  have,  from  Equations  (28a)  and  (33),  for 
a  control  line  lying  behind  the  lifting  line. 


The  corresponding  expression  for  w  is  therefore  taken  to  be 


^  i  \  g  cos  ;  -TT  ^  V'o  1  (35) 

00  00  ‘ 

D"0,2,4 

Note  that  these  series  are  unique  since  w(f^)  ■  w(-f^)  by  the 
choice  of  a  cosine  series,  so  that  w  is  specified  over  the  entire 
interval  -v  <  <  v  although  the  wing  span  covers  only  the 

interval  0  <  V'  ^  it  • 

Substitution  of  Equations  (34)  and  (35)  into  Equation  (32) 
yields  the  set  of  equations  to  be  solved  for  the  unknown  coefficients 
representing  the  circulation  distribution 

00 

7  -  ^  7^  sin  nf 

n-1,3,5 

along  the  lifting  line.  The  expanded  form  of  Equation  (32)  thus 
becomes 
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V4TrsV„^ 
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2  50  52  o  54 


X  idT/djlc 

4v-  V^s 
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(go  +  9 2  cos  2f^  +  cos  4^/^  +  ...)  +  sin  a 


(36) 


Now,  in  order  that  this  equation  be  satisfied  for  all  values  of 
V'q  it  must  be  satisfied  for  each  Fourier  component.  Therefore, 
if  we  let 


n-Y, 


'n  4Trsy_  sin  a 


and 


r*  .  (^r/<?£)g 


(37) 


4TrsV_  sin  a  2fAV  sin  a 
00  00 

the  following  set  of  equations  results : 

n  ■  1  n  -  3  n  -  5  n  -  7 
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These  equations  can  now  be  written  in  the  following  abbreviated 
form: 

00 

f.'y*“g.r*  +  l;  j“0 

nj'n  ^3 


where 


n** *  *  • 
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n* 1 , 3 , 5 


,  1 9  • 


f.“27r+5'a.;  j“0 

nj  2  n] 
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nD 


^  a 
2  ^nj 


4v  +  ^  a  . 
2  nj 


>  (39) 


f  .7  *  ■  g.r*  ;  j  *•  2,4,6,  ... 

n^fn  ^3  '  J  »  »  » 


2,4,6,  ... 


(40) 


The  numbers  denoted  by  f_ . 
which  can  be  used  to  obtain  the  loading  coefficients  7^^* 


are,  then,  universal  numbers 

for 

any  specified  separated  flow  pattern  as  given  by  the  coefficients 
g^ .  The  particular  type  of  flow  separation  assumed  here  is  shown 
in  Figure  1  and  the  corresponding  function  w(x,y)  is  given  by 
Equations  (10)  and  (12). 

In  the  asymptotic  formulation  for  large  distances  behind  the 
lifting  line,  all  the  coefficients  a^^  become  zero,  (see  Eqs. 


(25)  and 

(26))  and  the 

matrix 

of  the 

foregoing  infinite 

set 

equations 

yields 

the  following 

matrix 

of  coefficients: 
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Thus  it  is  clear  that  in  this  case  the  determinant  is  not  zero. 
Numerical  calculations  for  other  locations  behind  the  lifting  line 
indicate  that  this  is  in  fact  true  in  general.  Therefore,  there 
exists  no  nontrivial  solution  to  the  homogeneous  set  of  equations 
2  “  Oj  the  general  solution  of  Equations  (39)  for  the 

loading  coefficients  7^^*  can  be  represented  as 

V  "(v)  (vl 

o 

The  first  term  (V  * ^  is  the  solution  of  the  equations 

00 

E  ‘njV  -  1  <«> 

n*  1 , 3 , 5 

00 

E  13  -  <«) 

n- 1 , 3 , 5 

which  corresponds  to  no  flow  separation  (r*  ■  0) .  On  the  other 
hand,  solution  of  the  equations 

00 

Z!  ^nj'^n*  “  “  0>2,4,6...)  (45) 

n- 1 , 3 , 5 

which  corresponds  to  the  image  system  required  for  the  separated 
flow  (r*  ■  1).  Thus  for  a  given  value  of  f*,  the  general  solu¬ 
tion  (which  is  unique)  is  given  by  Equation  (42).  The  solution 
of  Equation  (45)  is  therefore  seen  to  be  the  basic  source  of  the 
change  in  the  aerodynamic  characteristics  of  the  wing  due  to 
separation. 
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4.2.5  Satisfying  the  Kutta  condition 

In  the  foregoing  solution,  the  value  of  r*  is  a  parameter 
which  has  not  been  specified,  and  a  solution  exists  for  each 
value  of  r*.  From  this  spectrum  of  solutions,  it  is  desired 
to  select  the  one  corresponding  to  the  value  of  r*  required 
on  the  basis  of  physical  considerations.  It  turns  out  that  it 
is  possible  to  make  the  determination  of  r*  on  the  basis  of  the 
Kutta  condition.  The  particular  condition  is  that  the  downwash 
just  off  the  side  edges  of  the  wing  be  finite  rather  than  infinite. 
Thus,  the  flow  at  the  side  edge  of  the  wing  seen  in  end  view  leaves 
the  side  edges  smoothly  and  (see  sketch)  does  not  turn  around  the 
sharp  edges  with  a  flow  singularity. 


The  mathematical  requirement  for  the  Kutta  condition  can  be 
established  by  examining  Equation  (22)  for  points  lying  outside 
the  span  of  the  vortex  sheet.  For  such  points,  by  making  use  of 
Equations  (24)  and  (25) ,  we  can  write  Equation  (22)  in  the  form 
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where  x  -  >  0,  since  we  need  to  satisfy  the  Kutta  condition 

behind  the  lifting  line  only.  Now,  integration  of  the  second  term 
yields,  for  points  outside  the  interval  of  integration  (Ref.  21) 


IT 

r  cos  Tit  dW  ^ 

r  , 

1  cos  V'  +  ^ 

^ 

for  1  ^ 

(47) 


Thus,  for  points  lying  outside  the  span  of  the  vortex  sheet  and 
behind  the  lifting  line,  the  upwash  w*/V^  becomes 


Hence,  the  singular  part  of  the  upwash  denoted  by  3'^®^ 

outside  the  span  is  given  by 


An  examination  of  this  result  reveals  several  very  interesting 
facts.  As  we  approach  the  side  edges  of  the  vortex  sheet  from  out¬ 
side  the  span  of  the  vortex  sheet,  the  downwash  for  each  Fourier 
component  has  a  square  root  singularity.  Furthermore  the  coeffi¬ 
cient  of  the  singularity  is  independent  of  x  so  that  the  singu¬ 
larities  have  the  same  strength  along  the  length  of  each  side  edge. 
It  follows,  therefore,  that  we  can  satisfy  the  Kutta  condition  along 
the  entire  length  of  each  side  edge  if 
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Z  "Yn  -  0  or  ^  7„*  -  0 

n- 1 , 3 , 5  n- 1 , 3 , 5 


(50) 


This  is,  therefore,  the  mathematical  formulation  of  the  Kutta 
condition^ . 

From  Equations  (42)  and  (50),  we  can  now  determine  the 
separated  vortex  strength  T*  which  is  necessary  to  satisfy  the 
Kutta  condition.  Thus,  substituting  Equation  (42)  into  Equation  (50), 
we  have 


or 
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00 

(v)„  *  '■*  I  (v),  -  0 


r* 


-  I 

— IgA-iPaJ 
00 
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n-1,3,5. 


(51) 


The  final  solution  of  the  complete  boundary-value  problem  of  a 
rectangular  wing  with  side-edge  separation  is  therefore  given  by 
Equations  (42)  and  (51),  so  that  the  loading  coefficients  of  the 
lifting  line  are 


^It  can  easily  be  shown  that  this  condition  is  precisely  the  same 
as  requiring  that  the  loading  7(1))  must  have  a  zero  slope  at 
the  wing  tips.  But  the  slope  of  the  loading  is  proportional  to 
the  strength  of  the  trailing  vortex  sheet,  so  the  Kutta  condition 
is  satisfied  along  the  side  edges  if  the  trailing  vortex  sheet 
has  zero  strength  at  the  edges. 
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■'n*  ■  ('^n*)„  -  (y«  X 


(52) 


It  is  clear  from  the  foregoing  analysis  that  there  is  only 
one  correct  solution  for  7^*  in  the  given  boundary-value  problem®. 
That  is,  the  solution  is  unique  (for  a  specified  0),  and  there  is 
only  one  correct  value  for  r*  which  satisfies  the  Kutta  condition. 
One  may  well  ask,  then,  what  is  the  meaning  of  the  so-called 
"unseparated"  solution  with  r*  ■  0?  Clearly,  this  solution  does 
not  satisfy  the  Kutta  condition  at  the  side  edges.  Mathematically, 
it  simply  represents  a  portion  of  the  complete  solution  of  the 
stated  boundary-value  problem.  But  physically  this  "unseparated" 
solution  actually  represents  the  solution  for  which  there  is  no 
side-edge  separation.  That  is,  the  side  edges  are  not  sharp,  and 
one  therefore  need  not  impose  the  Kutta  condition.  We  have  a 
different  boundary-value  problem,  then,  and  the  unseparated  solu¬ 
tion  therefore  corresponds  to  the  classical  wing  theory  of  Prandtl, 
which  applies  to  high  aspect  ratio  wings  at  low  angles  of  attack. 

4.2.6  Determination  of  the  shedding  angle  Q 

The  orientation  of  the  vortex  lines  which  comprise  a  free 
vortex  sheet  is  determined  by  the  average  of  the  velocity  vectors 
on  the  two  sides  of  the  sheet.  This  follows  from  the  Helmholtz 
law  that  a  free  vortex  must  follow  the  local  streamlines.  There¬ 
fore,  we  can  determine  the  inclination,  or  shedding  angle,  6,  of 
the  .ide-edge  vortex  sheets  by  calculating  the  total  velocity  com¬ 
ponents  just  above  the  wing  on  either  side  of  the  sheet.  Thus 
(see  sketch) , 

®Note  that  the  determination  of  each  loading  coefficient  for  the 
conplete  problem  requires  the  determination  of  all  of  the  loading 
coefficients  from  the  homogeneous  and  the  nonhomogeneous  solutions. 
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Since  w  and  u  are  both  functions  of  0,  the  determination  of 
6  will  require  an  iterative  process  in  which  a  value  of  6  is 
assumed,  the  right-hand  side  of  Equation  (53)  is  calculated  and 
compared  with  the  assumed  value  of  0,  the  right-hand  side  is 
recalculated  using  the  new  value  of  6,  and  so  on.  The  procedure 
is  analogous  to  that  employed  by  Bollay  (Ref.  1). 

Since  the  boundary  condition  of  no  flow  through  the  wing  has 
been  satisfied  only  at  one  chordwise  location  (i.e.,  at  the  control 
line),  we  shall  solve  Equation  (53)  for  9  at  the  control  line 
also.  At  this  location,  the  boundary  condition  requires  that 

^  +  y-  +  sin  o 
00  00 
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(54) 
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Now,  since  the  singular  part  of  w*  just  outside  the  span 
has  been  removed  by  imposing  the  Kutta  condition  of  Equation  (50), 
the  value  of  w*  at  y  ■  s  +  e  is  given  by  (see  Eq.  (28)) 


(55) 


and  it  is  noted  that  a  is  continuous  across  the  vortex  sheet. 

n 

The  value  of  w  just  beyond  the  wing  tip  can  be  found  simply  from 
Equation  (12c)  by  noting  that  as  we  pass  from  y“  s  -  e  to 
y  “  s  +  €  the  second  arc  tangent  in  that  equation  jumps  from 
0  to  2v,  while  the  first  arc  tangent  remains  at  ir.  That  is, 
the  last  term  of  Equation  (12c)  is  discontinuous  as  y  moves 
across  the  vortex  sheet,  and  we  find  that  Wj^  simply  changes  sign 
across  the  sheet.  Therefore,  we  can  write  w  just  outside  the 
span  in  the  form 


w  *  w 

y“s+e  y"s-€ 


-S-6 
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(56) 


Inasmuch  as  the  vortex  sheet  lying  in  the  plane  of  the  wing 
consists  of  a  single  vortex  line  ahead  of  the  control  line  and  a 
system  of  chordwise  trailing  vortices,  this  vortex  system  produces 
no  chordwise  velocity  component  at  the  control  line.  Therefore 

)  -(^)  ■“ 

y«S-€  y*S+€ 

Also,  it  can  be  seen  from  Appendix  A  that  the  chordwise  components 
of  velocity  due  to  the  side-edge  vortex  sheets  are  direqtly  related 
to  the  corresponding  upwash  in  the  following  manner: 

"r  ’  -^R  ® 


(58) 
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(Note  that  u  changes  sign  with  w  as  we  cross  the  vortex  sheet.) 
The  chordwise  velocity  component  produced  by  the  bound  portion  of 
the  separated  vortices  can  be  obtained  from  Appendix  A.  Thus, 
replacing  F  by  (dF/d^)  d|  and  integrating  over  the  chord,  we 
find  that 

c 

—  _  zs  dF  f  _ d£ _ 

ar  J  .  z»]  V(x  -  i)‘  *  4.=  +  z= 

Now,  integrating  this  expression  (Ref.  21)  and  taking  the  limit 
as  z  -*■  0  from  the  positive  direction,  to  find  the  induced  chord- 
wise  velocity  component  immediately  above  the  wing  surface,  we 
obtain 
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Finally,  substitution  of  Equations  (54)  through  (60)  into 
Equation  (53)  yields  the  following  expression  for  the  shedding 
angle: 
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(61) 


It  is  noted  that  Bollay  (Ref.  1)  obtained  a  value  of  0  equal 
to  half  the  angle  of  attack  for  the  limiting  case  of  vanishing 
aspect  ratio  at  small  angles  of  attack.  In  the  more  general  case 
(Eq.  (61)),  the  ratio  of  0/a  will  be  a  function  of  both  aspect 
ratio  and  angle  of  attack. 
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4.3  Normal  Force  and  Center  of  Pressure 

According  to  the  Kutta-Joukowski  law,  the  normal  force  acting 
on  the  lifting  line  of  circulation  strength  7(ri)  is  given  by 

+s 

Nj^  ■  pV^  cos  a  J  y  (t])  dt)  (62) 

-8 

so  that,  with  the  definitions  of  Equations  (13)  and  (14),  we  find 
that 


N,  ■  pV  cos 
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pV^  s  cos  a 
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■  2ir^pV^^s^  cos  a  sin  a  (j 

Thus,  only  7^*  contributes  to  the  normal  force,  whereas  the  higher 
harmonics  affect  only  the  span  load  distribution.  The  load  distri¬ 
bution  associated  with  7^*  is  elliptical,  but  its  magnitude  will 
be  affected  by  the  existence  of  the  higher  harmonics. 

For  the  unseparated  case  (F*  ■  0),  we  have  for  the  normal 
force  coefficient 
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(l/2)pV^2(2  sc) 


cos  a  sin  a 


(64) 


For  small  aspect  ratios,  the  lifting  line  and  the  line  of  co.*trol 
points  may  be  separated  by  many  semispans  so  that  (x  -  ^^)/s 
is  large.  Thus,  approaches  its  asyn^totic  value  of 

l/(2Tr)  and  we  have 

f  ^  ® 
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(65) 
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This  result  agrees  with  the  slender-body  result  for  the  normal 
force  coefficient  of  a  rectangular  wing  of  low  aspect  ratio  at 
low  angles  of  attack. 

The  total  normal  force  acting  on  the  wing  with  flow  separ¬ 
ation  is  the  sum  of  the  normal  force  on  the  lifting  line 
(including  the  image  system  (7  *)_)  and  the  normal  force  acting 
on  the  bound  vortices  of  the  separated  system  (which  have  uniform 
span  loading).  Thus,  using  Equations  (37),  (42),  (62),  and  (63), 
we  have  for  the  total  normal  force 
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The  total  normal  force  coefficient  is  therefore  given  by 


JL 


ipV„^(2sc)  2 


A  cos 


a  sin  +  r*^2Tr^7^*^  +  sj 


*  ^N_ 
u  S 


where 


C„  ■  o  A  cos 


a  sin  a  •  F*  (y ^ 


(67) 


(68) 


represents  the  normal  force  due  to  the  separated  vortices  and  their 
images . 

The  center  of  pressure  is  obtained  by  noting  that  the  load¬ 
ings  due  to  to  ('Yi*)5  both  act  on  the  lifting  line 

(i.e.,  at  the  quarter  chord),  whereas  that  due  to  F*  acts  at 
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the  half  chord,  since  dP/d^  was  taken  as  constant.  Thus  we  can 
write  at  once,  from  Equation  (67) 


(69) 


The  foregoing  solution  also  yields  a  precise  determination 
of  the  proper  location  of  the  control  line  if  the  present  theory 
is  to  yield  results  in  accordance  with  experiment  for  small  angles 
of  attack.  To  demonstrate  this  fact,  let  us  calculate  the  lift- 
curve  slope  at  zero  angle  of  attack  as  a  function  of  aspect  ratio 
and  position  of  the  control  line.  By  comparing  these  results  with 
those  from  experiment  or  exact  theory,  it  can  be  determined  where 
the  control  line  should  be  placed  as  a  function  of  aspect  ratio. 

The  lift-curve  slope  at  zero  angle  of  attack  has  been  computed 
as  a  function  of  aspect  ratio  using  Equation  (67).  Note  that 


0 


since  there  is  no  flow  separation  at  zero  angle  of  attack.  The 
results  are  shown  in  Figure  2  for  several  chordwise  positions  of 
the  control  line  with  the  lifting  line  placed  at  the  quarter- 
chord  position.  They  are  also  compared  with  the  mean  experimental 
curve  of  Scholz  (Ref.  5).  It  is  seen  that  locating  the  control 
line  at  the  three -quarter-chord  position  yields  excellent  agreement 
between  the  present  theory  and  experiment  for  unseparated  flows . 
-The  use  of  a  control  line  at  the  half-chord  position  yields  values 
of  lift-curve  slope  about  10  percent  low  at  an  aspect  ratio  of 
unity  and  25  percent  low  at  A  -  2.  These  results  are  in  substan¬ 
tial  agreement  with  the  classical  work  of  Weissinger  (Ref.  22) 
which  was  evidently  based  on  his  observation  that  the  influence  of 
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the  semi-infinite  trailing  vortices  varies  very  little  over  the 
airfoil  chord,  so  the  influence  of  the  bound  vortex  is  dominant 
in  determining  the  proper  placement  of  the  control  line.  The 
present  calculations  take  the  trailing  vortices  fully  into  account 
and  therefore  serve  to  substantiate  Weissinger's  conclusion. 

4.4  Dovmwash  at  the  Tail  Location 

In  order  to  investigate  longitudinal  stability  and  control 
characteristics  of  an  aircraft  configuration,  one  must  determine 
the  dcMiwash  at  the  location  of  the  horizontal  tail.  This  is  a 
relatively  straightforward  procedure  once  the  strengths  and 
positions  of  the  vortices  shed  from  the  wing  are  known.  Thus, 
the  upwash  at  any  point  in  space  is  given  by 

w(x,y,z)  ■  w*  +  w  +  sin  a 

and  the  expressions  for  w*  and  w  for  points  lying  in  the  plane 
of  the  wing  are  given  by  Equations  (12)  and  either  (28)  or  (48) , 
depending  upon  whether  the  point  in  question  lies  inside  or  out¬ 
side  of  the  wing  span. 

The  determination  of  the  upwash  w*  and  w  for  points  lying 
out  of  the  plane  of  the  wing  is  clearly  not  covered  by  Equations 
(12) ,  (28) ,  and  (48) ,  since  they  were  specialized  by  setting  z  ■  0. 
Therefore,  one  must  start  with  the  more  general  expressions  of 
Appendix  A  for  the  single  vortex  and  develop  expressions  analogous 
to  Equations  (12),  (28),  and  (48)  for  z  0.  The  procedure  is 
exactly  parallel  to  that  of  Section  4.2  of  the  present  report. 
However,  the  resulting  expressions  are  considerably  more  compli¬ 
cated  and  do  not  lend  themselves  to  the  same  method  of  evaluation.^ 


®Vidya  has  actually  carried  out  the  analysis  for  z  0,  but 
that  anedysis  has  not  been  checked  independently  and  has, 
therefore,  not  been  prograutied  for  ccxnputation.  The  analysis 
can,  however,  be  made  available  in  its  present  form. 
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4.5  Calculative  Procedure 

The  calculation  of  the  aerodynamic  characteristics  of  a  rec¬ 
tangular  wing  element,  using  the  analysis  developed  in  the  present 
report,  requires  the  solution  of  the  set  of  simultaneous  equations 
given  by  either  Equations  (38)  or  (39) .  This  solution  yields  the 
loading  coefficients,  strength  of  the  separated 

vortices,  r*,  which  in  turn  determine  the  normal  force  coefficient 
and  the  center  of  pressure. 

The  procedure  for  performing  the  numerical  calculations  is 
outlined  below: 

(a)  Select  the  number  of  terms,  n^,  to  be  used  in  the  series 
representation  of  the  loading  (i.e.,  the  nunJser  of  nonzero  harmonics 
to  be  carried  in  the  calculation) .  This  will  be  discussed  further 
in  Section  6. 

(b)  Calculate  the  Fourier  coefficients,  a^,  using  Equa¬ 
tions  (25)  and  (26)  with 

n  -  1,  3 ,  . . . ,  (2nj^  -  1) 

for 

0  <  f  o  <  ^ 

and 

X  -  ^ .  (x/c)  -  0.25 

^  m  ^ 

s  A/2 

where  x is  the  distance  of  the  control  line  behind  the  lead¬ 
ing  edge ,  expressed  in  chord  lengths . 

(c)  Calculate  the  Fourier  coefficients  a^^  where 

cos  (70) 

and 

j  ■  0,  2 .  (2nj^  -  2) 
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(d)  Determine  the  coefficient  matrix  of  the  set  of  simul¬ 
taneous  equations  (Eq.  (38)  or  (39))  using  Equation  (40). 

(e)  Assume  an  initial  value  of  the  shedding  angle,  Q,  for 
the  separated  vortex  system. 

(f)  Calculate  the  Fourier  coefficients,  g ^ ,  from 


(t)-i  / 


^  w(3V's,  ^q) 

■jTrr; - = -  cos  itO  di/ 

r*V  sin  a  -^^o 


(71) 


for 


j  *  0>  2^  •••>  (2 —  2 ) 


and 

X  ■  ^</’= 
s 

by  solving  Equations  (43)  and  (44). 
by  solving  Equation  (45) . 

(i)  Calculate  r*  using  Equation  (51) . 

(j)  Calculate  the  normal  force  coefficient,  C^^,  the  center 

of  pressure,  x/c,  and  the  new  shedding  angle,  0,  using  Equa¬ 
tions  (67),  (69),  and  (61),  respectively. 

(k)  Using  this  new  value  of  0,  repeat  steps  (f)  through  (j) 
and  compare  the  two  values  of  C^j.  If  the  desired  accuracy  in 

has  not  been  obtained,  return  to  step  (f)  and  repeat  until  it  is 
obtained . 


(g)  Determine 

(h)  Determine  fy 

V"  A 
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5.  AERODYNAMIC  THEORY  OF  SWEPTBACK  WINGS  WITH  LEADING-EDGE 
SEPARATION 

5.1  Construction  of  Mathematical  Model 

Having  developed  a  theoreticarl  analysis  for  the  rectangular 
wing  element  with  side-edge  flow  separation  which  satisfies  the 
proper  conditions  of  no  flow  through  the  wing  and  finite  veloc¬ 
ities  at  the  side  edges,  we  can  now  proceed  to  apply  the  theory 
to  a  wider  variety  of  planforms  in  the  following  manner.  Since 
the  Kutta  condition  of  finite  velocity  at  the  side  edges  is  actu¬ 
ally  satisfied  for  all  values  of  x  (see  Eq.  (50)),  the  only 
singularities  in  velocity  in  the  foregoing  theory  are  along  the 
lifting  line  itself.  Therefore,  o^ther  planforms  can  be  constructed 
by  using  rectangular  wings  of  various  aspect  ratios  as  elementary 
building  blocks,  provided  that  the  trailing  edge  is  straight.  This 
representation  will  enable  us  to  treat  rather  general  leading-edge 
shapes,  as  shown  in  the  sketch  below. 


lifting  line  (at  1/4 
chord)  of  rectangular 
wing  element 
control  line  of  rectan¬ 
gular  wing  element 
(at  3/4  chord) 


We  shall  now  place  a  lifting  line  at  the  quarter  chord  of  each 
rectangular  wing  element  and  proceed  to  satisfy  the  boundary  con¬ 
dition  of  no  flow  through  the  wing  surface  at  the  corresponding 
3/4-chord  lines,  accounting  for  the  downwash  from  all  of  the  wing 
elements  and  their  associated  shed  vortices.  This  will  require 
that  we  be  able  to  calculate  the  downwash  produced  along  the 
3/4  chord  of  any  rectangular  wing  element  by  the  vortex  system 
associated  with  any  other  wing  element.  It  is  clear,  then,  from 
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the  foregoing  sketch  that  for  a  sweptback  wing  we  must  be  able  to 
calculate  the  downwash  anywhere  in  the  plane  of  the  wing  as  pro¬ 
duced  by  any  of  the  rectangular  wing  elements.  In  particular, 
we  shall  require  the  downwash  outside  the  span  of  the  lifting  line 
itself.  For  this  purpose,  we  shall  have  recourse  to  Equation  (48) 
which  was  used  in  satisfying  the  Kutta  condition. 

Consider  now  the  case  of  a  sweptback  wing  which  is  to  be  made 
up  of  n^  rectangular  wing  elemints  of  equal  chord  c^^  ■  °c/”w 
and  whose  semispans  are  selected  such  that  the  area  of  each  rectan¬ 
gular  element  is  equal  to  the  area  it  replaces  (see  sketch) 


c . 
1 


Note  that  the  distances  x.  and  will  be  measured  rearward 

1  th  " 

from  the  leading  edge  of  the  k  wing  element. 

The  boundary  condition  to  be  satisfied  on  each  of  the  n^ 

rectangular  wing  elements  is  that  the  normal  velocity,  (w) . . , 

th 

due  to  the  separation  vortices  of  the  wing  itself,  the  i  wing, 
plus  the  normal  velocity,  (w*)..,  due  to  the  lifting  line  of 
the  i  wing  and  its  trailing  sheet  must  exactly  cancel  the 
component  of  the  free  stresun  normal  to  the  wing  (V^  sin  a)  plus 
the  sum  of  the  normal  components  induced  by  the  vortex  systems 
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o£  all  of  the  other  wings , 


(w)i]^. 


where 


**_  _  ** 

w  r 

k-1  k-i 


'">ik 


The  boundary  condition  can,  therefore,  be  written  in  the  form 

11  11  k*i  L  ik  ik 

where  this  condition  must  be  satisfied  for  all  values  of 


i(i  ■  1,  2,  ...,  n^) .  Thus,  a  set  of  n^  simultaneous  equa¬ 
tions  is  generated.  Now  if  the  loading  coefficients  and  the 
separated  vortex  strengths  on  all  of  the  rectangular  wing  ele¬ 
ments  except  the  i^^  wing  are  known ,  and  if  the  summation  on 
the  right-hand  side  of  Equation  (73)  can  be  replaced  by  an  even 
cosine  series ,  then  the  boundary  value  problem  can  be  solved  in 
the  same  manner  as  is  presented  in  Section  4.2.4  for  the  rectan¬ 
gular  wing  alone. 

Initially,  in  the  solution  of  the  sweptback  wing  problem, 
the  vortex  strengths  and  the  loading  coefficients  are  not  known 
for  any  of  the  n  wings.  However,  an  initial  guess  can  be  made 
by  treating  each  wing  as  an  isolated  rectangular  wing  and  using 
the  method  of  Section  4  to  obtain  initial  values  for  the  loading 
coefficients,  (7^*^  separated  vortex  strengths, 

for  each  of  the  n^  ^wings.  Then  an  iterative  technique  can  be 
used  to  solve  the  set  of  n^  equations.  Since  the  sweptback 
wings  being  treated  in  the  present  analysis  are  restricted  to 
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those  which  are  synunetrical  about  the  midspan,  the  midspans  of 
the  rectangular  elements  are  aligned  and  thus  the  downwash  pro¬ 
duced  on  any  element  by  any  other  elements  is  symmetrical  about 
the  midspan  of  the  i^^  element.  Therefore,  the  summation  can 
be  represented  by  an  even  cosine  series  of  the  form 


f — '  ik  ik  ._*■ — 1 


k=i 

k/^i 


(74) 


D“0,2,4 


where 


~  IT  £.  ^  TT 


Now,  by  a  development  similar  to  that  of  Section  4.2.4  the 
set  of  simultaneous  equations  to  be  solved  may  be  written  as 


w 


n“ 1 , 3 , 5 
00 


)  (75) 


X  ('nj),,(V).  -(^Ou  "i*  ‘ 


2,4,... 


n* 1 , 3 , 5 


There  are  n^  sets  of  these  simultaneous  equations  to  be  solved 
by  the  iterative  technique.  Equations  (75)  are  the  counterpart 
of  Equations  (39)  for  the  rectangular  wing  and  may  be  solved  in 
the  same  manner  by  splitting  the  solution  for  y^*  into  an 
unseparated  solution  (r^*"0)  and  a  solution  due  to  separation, 
as  described  in  Section  4.2.4. 

In  order  to  obtain  the  coefficients  b^ ,  we  shall  have  to  cal 
culate  the  upwash  w^j^  produced  at  the  3/4  chord  of  the  i^li  wing 
element  by  the  vortex  system  associated  with  the  k^^  wing  element 
For  this  purpose,  we  shall  need  the  distances  from  the  ktb  lift¬ 
ing  line  to  the  i^^  control  line,  given  by 


53 


(76) 


for  the  determination  of  w..  *,  and  the  distance  from  the  lead- 

til  til 

ing  edge  of  the  k  wing  element  to  the  i  control  line, 
given  by 


(77) 


for  the  determination  of  w..  .  Note  that  both  of  these  distances 

xk 

are  expressed  in  semispans  of  the  wing  element  whose  vortex  system 
produces  the  downwash  (i.e.,  of  the  k^^  wing  element). 

Now  from  Equation  (28)  we  can  write  for  points  lying  within 
the  span  of  the  lifting  line  (i.e.,  for  -1^ 


_  ^j-^k 

V  sin  a  V  Si  ’ 


w.,  * 
ik 


..y  L .)  .  I  a 

„L^  "4  ♦'o  2  ”k\  “k  °/ 

n*r, 3  L  J 


for 


Xi  -  ^k  >  0 


(78) 


■  ^  2  XV  ®k 

n*i , 3  L 

for  x^  -  <  0 


(79) 


where 


lb  ■  COS 

^  o 


■■  (•  t) 
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Por  points  lying  outside  the  span  of  the  lifting  line  (i.e. ,  for 
1  Equation  (28)  is  modified  in  accordance  with 

Equation  (48)  to  yield 


W  ^ 

''ik 


V  sin  a 
00 


n 


for  Xi  -  <  0 


(81) 


Similarly,  we  obtain  the  required  expression  for  w..  (the  up- 

til 

wash  produced  at  the  i  control  line  by  the  separated  vortex 
system  associated  with  the  k^^  wing  element)  by  replacing  x 
by  x^,  s  by  Sj^,  and  c  by  c^  in  Equation  (12).  With  the 
resulting  equation  and  Equations  (78)  through  (81)  ,  we  are  now  in 
a  position  to  set  up  an  iterative  procedure  for  determining  the 
loading  on  the  sweptback  wing  with  leading-edge  separation.  The 
detailed  procedure  will  be  outlined  in  Section  5.5. 
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It  should  be  noted  that  the  same  method  may  be  applied  to  a 
low  or  intermediate  aspect  ratio  rectangular  wing  by  representing 
it  as  a  chordwise  distribution  of  higher  aspect  ratio  rectangular 
wing  elements.  The  purpose  of  such  a  procedure  in  these  cases 
would  be  to  extract  more  detailed  information  regarding  chordwise 
load  distribution  and  center  of  pressure  than  can  be  obtained  by 
using  a  single  wing  element. 

5.2  Determination  of  Shedding  Angle,  Q.,  on  i^^  Rectangular 
Wing  Element  ^ 

The  determination  of  the  angle,  6^,  at  which  the  separated 
vortices  are  shed  from  each  of  the  rectangular  wing  elements  is 
similar  to  that  described  in  Section  4.2.6  for  the  rectangular 
wing  alone.  However,  in  addition  to  the  velocities  induced  by  the 
ith  wing  itself,  we  must  now  consider  the  velocities  induced  by 
all  of  the  other  wing  elements  making  up  the  swept  wing. 

The  general  expression  for  the  shedding  angle  on  the  i^^ 
wing  is  given  by  (Eq.  (53)) 


tan  6 . 
1 


w. 


av 


av 


2  (^y-s^-e  ^y-Sj^+e) 
2  ('^ySi-e  “y-s^+e) 


(82) 
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This  equation  is  solved  by  an  iterative  procedure,  as  was  Equation  (61)  for  the 
single  rectangular  wing  element. 
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5.3  Normal  Force,  Center  of  Pressure,  and  Span  Loading 

The  total  normal  force  on  the  sweptback  wing  is  sinply  the 
sum  of  the  normal  forces  on  all  of  the  wing  elements,  so  that  the 
total  normal  force  coefficient  is,  from  Equations  (66)  and  (67), 


w 


E  -i 


JL 


1  2 


1»1 


pv.-s 


00  O 


W 


2  pv."  E 


1-1 


"^W  2 

47rn^  sin  a  cos  a  ^  ^  C%)[2  GxO.  ^ 


"w 


(87) 


E  ^ 


The  center  of  pressure  is  found  by  summing  the  pitching  moment 
about  the  wing  leading  edge  and  dividing  by  the  total  normal  force. 
That  is , 


X 


E  "A 

i»i _ 


i-i 


(88) 


Thus,  noting  that  the  loading  due  to  (7  *)  acts  at  the  quarter 
chord  of  the  i^^  wing  element  whereas  r^*  acts  at  the  midchord 
of  the  element,  we  have 
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w  c  r 

Z  N.3f.  -  4ir  ^  pvj  sin  a  cos  a  ^  s  |  (y/),  (i  -  J) 
i“  1  i“  1  ^  ^ 

2r/  (i  -  0]  (89) 

Therefore,  the  center  of  pressure  for  the  sweptback  wing  is  given 
by 


X  _  l-l 

c 


n 

w  xg.  vf  r  1 

n  yf—  )  +2r.* 

w  ^  \^3  J  1_2  1  J 


(90) 


The  span  loading  (in,  say,  pounds  per  foot)  is  obtained  by 
expressing  the  quantity  c^c  at  each  spanwise  station,  where  c^^ 
is  the  section  lift  coefficient.  Equating  two  expressions  for 
the  lift  on  an  elementary  chordwise  strip  (see  sketch) , 


we  have 


c^c  dy  •  i 


pV„r  dy 


(91) 
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where  T  is  the  total  circulation  on  the  chordwise  strip.  Thus 


c^c 


-  2L 


(92) 


In  the  present  mathematical  model,  the  total  circulation  round 
any  chordwise  strip  is  given  by 


r(y) 


(93) 


where  the  summation  is  carried  over  the  wing  elements  whose  semi¬ 
span  is  greater  than  the  value  of  y  in  question.  Thus,  since 
7(y)  in  the  present  analysis  is  represented  by  a  Fourier  sine 
series.  Equations  (92)  and  (93)  combine  to  yield  the  following 
expression  for  the  non-dimensional  span  loading  on  a  swept  wing. 


=i=(y)  _ 

O  no  O 


W 


— - — *•  sin  nf. 
n 


n»i ,  3 


(94) 

where  ■  cos"^(-y/s^)  and  the  summation  over  i  is  carried 
over  the  wing  elements  whose  semispan  is  greater  than  y. 
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5.4  Downwash  at  the  Tail  Location 

The  calculation  of  the  dovmwash  behind  a  swept  wing  is  again 
a  matter  of  collecting  the  appropriate  expressions  for  w*  and 
w.  That  is ,  the  upwash  at  any  point  is  again  given  by 

w(x,  y,  z)  -  w*  +  w  +  sin  a 

just  as  in  the  rectangular  wing  case,  but  here  we  must  account  for 
the  contribution  of  every  wing  element  comprising  the  swept  wing. 
Thus ,  we  have 

"w 

“  Z  '"ik* 

k-i 

and 

'  ■  I 

k-i 

where  the  distance  Xj^  appearing  in  the  expressions  for 
and  w^^  is  now  the  distance  from  the  leading  edge  of  the  k^^ 
wing  element  to  the  point  in  question.  For  points  lying  in  the 
plcuie  of  the  wing,  the  expressions  for  are  given  by  Equa¬ 

tions  (78)  to  (81) ,  depending  upon  whether  the  point  in  question 
lies  inside  or  outside  of  the  span  of  the  k^^  wing  element. 

The  corresponding  expressions  for  w^^  are  given  by  Equations  (12) 
if  we  replace  x  by  x^,  s  by  s^^,  and  c  by  Cj^. 

For  points  lying  out  of  the  plane  of  the  wing,  see  Footnote  9, 
page  46. 


62- 


5.5  Calculative  Procedure 

The  detailed  procedure  for  performing  the  numerical  calcula¬ 
tions  is  outlined  below: 

(a)  Select  the  number  of  wing  elements,  n^,  and  the  number 
of  harmonics,  n^,  to  be  carried  in  the  calculation.  This  point 
will  be  discussed  in  Section  6.2  where  convergence  with  respect  to 
these  parameters  is  investigated. 

(b)  Assume  an  initial  value  of  the  shedding  angle, 
e. (i  =  1,  2  ...  n  ),  for  each  wing  element. 

(c)  Calculate  the  coefficient  matrix  for  each  wing, 

treating  each  as  an  isolated  wing,  using  the  method  outlined  in 
steps  (b)  through  (d)  of  Section  4.5. 

(d)  Calculate  the  Fourier  coefficients  (g.;)..  as  outlined 

li 

in  step  (f)  of  Section  4.5. 

(e)  Obtain  initial  values  of  the  loading  coefficients, 

(y  *)  and  (7  *)  ,  for  each  wing  element  by  treating  each 

u  i  "s  i 

wing  as  an  isolated  rectangular  wing  and  solving  Equations  (43) , 
(44),  and  (45). 

(f )  Calculate  for  each  isolated  wing  by  Equation  (51) . 

(g)  Calculate  for  i  -  1;  k  -  2,  ...,  n^  using 

Equation  (76)  and  Equations  (78),  (79),  (80),  or  (81)  as  appro¬ 
priate  with  the  values  of  from  step  (e)  . 

(h)  Calculate  Wj^j^  for  i  ■  1;  k  ■  2,  . .  . ,  n^  using  Equa¬ 
tions  (12)  and  replacing  x  by  x^  (Eq.  (77)),  s  by  Sj^,  and 
c  by  Cj^. 

(i)  Determine  coefficients  (74), 

using  the  results  of  steps  (g)  and  (h) . 

(j)  Solve  Equations  (75)  for  the  unseparated  loading  coef¬ 
ficients  (Yn„*)  .  by  setting  r.*  •  0  for  i  ■  1. 

u  ^  X 

(k)  Solve  Equations  (75)  for  the  loading  coefficients  due 

to  separation  ('Yno*^  ■  ^  setting  Fj*  “  1  and  omitting  the 

remainder  of  the  right-hand  side  of  Equations  (75)  for  i  >  1. 
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(l)  Replace  results  of  step  (e)  with  results  of  steps  (j) 
and  (k) . 

(m)  Determine  r^*  for  i  -  1  by  using  Equation  (51) . 

This  completes  the  first  iteration  on  (7n*)  • 

i 

(n)  Replace  results  of  step  (f)  with  results  of  step  (m) . 

(o)  Repeat  steps  (g)  through  (n)  for  i  ■  2,  3,  n^ 

with  i  k  in  steps  (g)  and  (h) .  This  completes  the  first 
iteration  for  all  wing  elements. 

(p)  Calculate  normal  force  coefficient  and  center  of  pres¬ 
sure  using  Equations  (87)  and  (90). 

(q)  Repeat  steps  (g)  through  (p)  until  desired  accuracy 
on  normal  force  is  obtained. 

(r)  Calculate  0. (i  ■  1,  2,  ...,  n  )  from  Equation  (86). 

X  w 

(s)  Repeat  steps  (d)  through  (r)  until  the  desired  accuracy 
on  6^  is  obtained. 
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6.  NUMERICAL  CALCULATIONS  AND  COMPARISON  WITH  EXPERIMENT 

This  section  of  the  report  will  describe  the  results  of  the 
calculations  wkich  were  made  using  the  analysis  presented  in  Sec¬ 
tion  4  for  the  rectangular  wing  and  that  presented  in  Section  5 
for  the  swept  wing.  The  rectangular  wing  calculations  will  be  dis¬ 
cussed  first,  since  this  wing  element  is  the  building  block  for  the 
swept  wing  analysis .  Convergence  of  the  rectangular  wing  solution 
is  therefore  necessary  if  the  swept  wing  solution  is  to  converge. 

The  rectangular  wing  calculations  were  carried  out  on  an  IBM  1620 
digital  computer,  whereas  the  swept  wing  calculations  were  performed 
on  an  IBM  7094  computer. 

6.1  Rectangular  Wing 

Two  investigations  were  made  using  the  rectangular  wing  analy¬ 
sis  presented  in  Section  4.  The  first  of  these  was  to  determine 
the  number  of  terms  (nonzero  Fourier  harmonics)  required  in  the 
series  representation  of  the  loading  in  order  to  insure  convergence 
on  values  of  the  normal  force  and  center  of  pressure.  The  second 
investigation  was  to  determine  the  angle  B  at  which  the  separated 
vortices  must  be  shed  in  order  for  the  vortices  to  follow  the  local 
streamlines  in  the  vertical  plane  containing  the  wing  tip.  After 
these  studies  were  conpleted,  systematic  calculations  of  the  normal- 
force  coefficient  and  center  of  pressure  were  carried  out  for  a 
range  of  aspect  ratios  and  angles  of  attack,  and  the  results  have 
been  compared  with  available  experiment.  In  all  of  these  calcula¬ 
tions,  the  rectangular  wing  is  represented  by  a  single  wing  element, 
since  we  wish  to  investigate  the  convergence  of  the  basic  building 
block  of  the  theory. 

6.1.1  Convergence  with  number  of  harmonics 

Table  II  presents  the  results  of  rectangular  wing  calculations 
in  which  the  number  of  nonzero  harmonics  in  the  loading  was 
varied  from  2  to  14,  for  two  angles  of  attack,  over  a  range  of 
aspect  ratios.  For  these  calculations,  the  control  line  was  fixed 
at  the  3/4-chord  location  and  the  shedding  angle  B  was  held  fixed 
at  half  the  angle  of  attack.  It  can  be  seen  that  the  values  of 
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C^,  x/c,  and  F*  converge  rapidly  and  that  8  nonzero  harmonica 
evidently  suffice  for  1 -percent  accuracy  in  even  at  an  aspect 

ratio  of  20.  The  values  of  }^c  and  F*  converge  somewhat  more 
slowly,  as  seen  from  Tables  11(b)  and  II  (c),  and  it  can  be  seen 
that  the  rate  of  convergence  decreases  as  the  aspect  ratio  is 
increased.  (This  point  is  significant  in  regard  to  the  swept  wing 
solution.)  Also,  it  can  be  seen  that  convergence  is  more  rapid 
at  the  higher  angles  of  attack,  and  it  is  evident  that  the  solution 
will  have  convergence  difficulties  at  vanishing  angles  of  attack, 
since  the  shed  vortices  tend  to  coalesce  into  the  wing  tip. 

The  significance  of  the  increase  in  the  number  of  harmonics 
required  with  increasing  aspect  ratio  is  illustrated  by  Figure  3, 
which  can  be  used  as  a  guide  for  the  swept  wing  solution.  For  the 
rectangular  wing  itself,  this  feature  of  the  solution  is  of  little 
significance,  since  one  would  logically  return  to  linear  theory  at 
very  high  aspect  ratios.  However,  for  a  swept  wing,  even  of  low 
aspect  ratio,  the  aspect  ratio  of  the  trailing  element  increases 
as  the  number  of  wing  elements  is  increased.  For  example,  if  a 
delta  wing  of  aspect  ratio  3  is  represented  by  10  rectangular 
wing  elements,  the  trailing  element  has  an  aspect  ratio  of  14.25. 
Thus,  Figure  3  may  be  used  as  a  guide  for  estimating  the  number 
of  harmonics  one  might  need  to  use  in  the  swept  wing  solution.  It 
must  be  borne  in  mind,  however,  that  the  mutual  interference  among 
wing  elements  may  require  still  larger  numbers  of  nonzero  harmonics 
for  convergence  of  the  swept  wing  solution. 

6.1.2  Detailed  loading  coefficients  for  a  specific  case 

The  detailed  convergence  of  the  mathematical  solution  is 
perhaps  best  illustrated  by  listing  all  of  the  loading  coefficients 
for  a  particular  rectangular  wing  element.  For  this  purpose,  we 
shall  select  a  rectangular  wing  of  aspect  ratio  3.5  at  20°  angle 
of  attack,  with  the  shedding  angle  6  fixed  at  half  the  angle  of 
attack  and  the  control  line  located  at  the  3/4~chord  position. 

The  reason  for  selecting  this  particular  case  is  that  it  represents 
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the  fourth  wing  element  of  a  delta  wing  of  aspect  ratio  2 .  We 
shall,  therefore,  be  able  to  compare  the  loading  coefficients 
on  the  isolated  wing  element  with  those  in  which  the  interference 
from  other  elements  is  included. 

The  loading  coefficients  for  the  selected  case  are  listed 
in  Ted>le  III ,  and  it  is  noted  that  19  nonzero  harmonics  have  been 
carried  in  the  calculation,  to  correspond  with  the  delta  wing 
calculations  of  Section  6.2.4.  It  can  be  seen  from  Table  III 
that  the  loading  coefficients  associated  with  the  higher  harmonics 
(n  large)  are  several  orders  of  magnitude  smaller  than  the  leading 
coefficients  so  that  one  need  not  have  carried  19  harmonics  for 
this  case.  In  fact,  as  was  pointed  out  in  the  previous  section 
(Fig.  3),  6  harmonics  was  a  sufficient  number  foi^  the  selected 
case . 

The  separated  vortex  strength  r*  and  the  normal  force  and 
center  of  pressure  corresponding  to  the  loading  coefficients  of 
Table  III  are  obtained  by  using  Equations  (51) ,  (67) ,  and  (69) . 

The  resulting  values  are 

r*  -  0.06307 

C„  -  1.4450 
N  . 

3^c  -  0.4043 

6.1.3  Shedding  angle  calculations 

Calculations  have  been  carried  out  to  determine  the  shedding 
angle  0  of  the  separation  vortices  by  the  iterative  procedure 
of  Section  4.2.6,  and  the  results  are  plotted  in  Figure  4.  It  can- 
be  seen  that  the  shedding  angle  ratio  d/o  increases  with  aspect 
ratio  and  decreases  with  angle  of  attack.  Furthermore,  it  is 
noted  that  the  shed  vortices  are  found  to  lie  above  the  free-stream 


direction  for  aspect  ratios  above  2  for  angles  of  attack  below 


-67- 


15°.  This  result  is  known  to  be  unrealistic  because  of  the  down- 
wash  which  must  be  produced  by  a  lifting  wing  and  which  must  cause 
the  trailing-edge  vortices  to  drift  below  the  free-stream  direction. 
On  the  other  hand,  it  must  be  borne  in  mind  that  the  assumed  mathe¬ 
matical  model  does  not  permit  rolling  up  of  the  side-edge  vortex 
sheets.  This  process  would  naturally  tend  to  drive  the  shed  vortices 
downward  and  result  in  a  smaller  shedding  angle  of  the  rolled-up 
vortex  cores. 

6.1.4  Effect  of  control  line  location 

Because  of  the  above  result,  and  because  the  placement  of  the 
control  line  (where  0  is  calculated)  at  the  3/4-chord  location 
was  based  on  matching  the  lift-curve  slope  at  0=0,  it  was  con¬ 
sidered  worthwhile  to  investigate  the  effect  of  control  line  loca¬ 
tion  on  the  calculated  shedding  angle  6.  Therefore,  0  was 
recomputed  with  the  control  line  placed  at  the  65  and  85  percent 
chord  locations ,  and  the  results  are  presented  in  Figure  5 .  It  can 
be  seen  that  the  effect  of  control  line  location  on  the  shedding 
angle  was  rather  small,  and  that  the  new  control  line  locations 
failed  to  yield  more  realistic  values  of  0.  Therefore,  it  was 
concluded  that  the  control  line  should  remain  at  the  3/4-chord 
location  ,  in  order  to  insure  the  correct  lift-curve  slope  at  small 
angles  of  attack,  but  that  the  values  of  the  shedding  angle  0 
determined  by  the  iterative  technique  of  the  present  analysis 
(using  a  single-wing  element)  are  evidently  unrealistic.  This 
fact  is  most  likely  associated  with  the  use  of  a  lifting  line 
rather  than  a  lifting  surface  and  with  the  fact  that  the  shedding 
angle  is  determined  only  at  the  location  of  the  control  line  at 
which  the  boundary  condition  is  satisfied.  This  situation  might 
be  expected  to  improve  if  one  were  to  use  a  number  of  wing  elements 
of  higher  aspect  ratio  to  represent  the  rectangular  wing.  This  has 
not  been  done . 
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6.1.5  Normal  force  and  center  of  pressure 

Systematic  calculations  have  been  preformed  to  determine  the 
variations  of  normal  force  and  center  of  pressure  on  rectangular 
wings  at  various  angles  of  attack  for  a  wide  range  of  aspect  ratios 
(using  a  single  wing  element) ,  and  these  results  are  presented  in 
Figure  6.  For  these  calculations,  the  control  line  is  at  the 
3/4-chord  location,  and  the  value  of  d  is  that  determined  by 
iteration.  For  a  given  angle  of  attack,  the  normal-force  coeffi¬ 
cient  is  found  to  increase  with  increasing  aspect  ratio  while  the 
center  of  pressure  moves  forward.  For  a  given  aspect  ratio,  the 
normal-force  coefficient  increases  nonlinearly  with  angle  of  attack, 
and  the  center  of  pressure  is  seen  to  be  relatively  insensitive  to 
angle  of  attack. 

6.1.6  Comparisons  with  experiment 

The  calculated  variation  of  normal  force  with  angle  of  attack 
(for  a  single  rectangular  wing  element),  is  compared  with  experi¬ 
mental  data  in  Figure  7(a)  for  rectangular  wings  of  aspect  ratios 
from  0.2  to  5.0.  It  can  be  seen  that  the  theoretical  normal  force 
agrees  with  experiment  within  about  10  percent  over  the  angle-of- 
attack  range  shown.  A  notable  exception  is  the  aspect  ratio  5 
wing  beyond  the  stall,  and  it  is  recalled  that  the  present  theoret¬ 
ical  model  does  not  account  for  leading-edge  stall  on  rectangular 
wings.  In  a  broader  sense,  it  may  be  said  that  the  theory  fails 
if  the  flow  becomes  unsteady.  Thus,  for  all  aspect  ratios,  the 
theory  will  fail  above  some  angle  of  attack,  since  all  flat  plates 
at  90°  angle  of  attack  produce  an  unsteady  wake,  whether  the  vortices 
in  question  arise  from  leading,  trailing,  or  side  edges.  However, 
the  angle  of  attack  for  the  onset  of  unsteady  shedding  will  be  lower 
for  the  higher  aspect  ratios,  because  of  the  orientation  of  the 
shed  vortices  (spanwise  rather  than  chordwise).  That  is  to  say, 
vortices  which  are  aligned  nearly  normal  to  the  free  stream  produce 
an  unsteady  flow,  while  those  aligned  nearly  with  the  free  stream 
produce  a  steady  flow.  Thus,  as  seen  in  Figure  7(a),  the  angle-of- 
attack  range  over  which  the  theory  is  valid  will  be  larger  for  the 
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lower  aspect  ratios.  In  Reference  1,  Bollay  showed  that  his  theory 
agreed  with  experiment  for  angles  of  attack  up  to  45°  for  an  aspect 
ratio  of  1/30.  The  present  theory  is  in  close  agreement  with  Bollay* s 
theory  for  vanishing  aspect  ratios,  at  least  insofar  as  normal  force 
is  concerned,  despite  the  assumptions  of  a  single  lifting  line  and 
a  constant  separated  vortex  strength  along  the  chord. 

In  Figure  7(b),  the  calculated  center  of  pressure  is  presented 
for  a  variety  of  rectangular  wing  elements.  It  is  seen  that  the 
present  theory  predicts  centers  6f  pressure  which  are  rather  far 
rearward  (toward  the  midchord  location).  Experimentally  (Ref.  4), 
it  is  found  that  the  center  of  pressure  remains  closer  to  the  quarter 
chord  for  rectangular  wings  of  all  aspect  ratios .  It  is  interesting 
to  note  that  Bollay 's  theory  (Ref.  1)  places  the  center  of  pressure 
at  the  quarter  chord  for  all  aspect  ratios  by  assuming  the  chordwise 
load  distribution  to  be  that  of  a  wing  of  infinite  aspect  ratio. 

The  present  theory,  on  the  other  hand,  places  the  lifting  line  at 
the  quarter  chord  but  places  the  center  of  pressure  of  the  separated 
vortices  at  the  midchord  by  assuming  the  separated  vortices  to  be 
uniformly  distributed  in  the  chordwise  direction.  This  restriction 
can  be  removed  by  employing  a  number  of  rectangular  wing  elements 
to  represent  the  rectangular  wing  and  then  solving  for  the  chordwise 
distribution  of  loading  and  shed  vorticity. 

For  all  of  the  above  results  using  a  single  rectangular  wing 
element,  the  shedding  angle  of  the  separated  vortices  has  been  cal> 
culated  by  the  iterative  procedure  described  in  Section  4.2.6.  It 
will  be  recalled  (see  Fig.  4)  that  these  calculated  shedding  angles 
appear  to  be  unrealistically  high,  particularly  at  high  aspect 
ratios  and  at  low  angles  of  attack.  Nevertheless,  the  agreement 
with  experimental  normal  forces  on  rectangular  wings  is  quite  satis¬ 
factory  even  in  these  ranges . 

Comparisons  with  experimental  values  of  the  shedding  angle 
are  of  questionable  validity,  since  the  mathematical  model  is  clearly 
unrealistic  in  this  detail.  That  is  to  say,  the  vortex  sheet  shed 
from  the  side  edge  of  a  rectangular  wing  actually  rolls  up  into  a 
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single  vortex  core  which  lies  above  the  wing  surface.  It  is  the 
position  of  this  core  which  is  measured  experimentally.  The  only 
known  measurements  of  this  kind  (contained  in  Ref.  23)  indicate  a 
considerable  (nonlinear)  variation  in  the  orientation  of  this  vortex 
core,  both  with  chordwise  location  and  angle  of  attack  for  a  rec¬ 
tangular  wing  of  aspect  ratio  2.  In  any  case,  it  would  not  be 
reasonable  to  expect  that  the  theoretical  shedding  angle  at  the 
side  edge  which  renders  the  mathematical  solution  of  the  problem 
unique  for  the  selected  theoretical  model  (with  no  rolling  up) 
should  agree  with  experimental  measurements  of  the  orientation  of 
the  rolled-up  vortex  core.  In  fact,  it  is  clear  that  rolling  up 
of  the  side-edge  vortex  sheets  would  cause  the  uppermost  shed 
vortices  to  be  depressed  toward  the  wing  surface.  Therefore,  the 
experimental  angle  6^  of  the  vortex  core  relative  to  the  plate 
should  logically  be  smaller  than  the  theoretical  angle  Q  of  the 
side-edge  sheets.  That  is 

6  >  e 

^  c 

It  is  of  some  interest  to  note  that  in  Bollay's  theory  (Ref.  1) , 
the  theoretical  values  of  shedding  angle  may  well  be  more  realistic, 
since  he  satisfied  the  boundary  condition  (of  no  flow  through  the 
wing)  in  the  mean  rather  than  at  a  specific  chordwise  location,  and 
also  since  he  assumed  a  more  realistic  chordwise  variation  of  shed 
vorticity  than  the  uniform  distribution  of  the  present  theory  for  a 
single-wing  element.  Bollay's  approach  yielded  shedding  angles  not 
far  from  half  the  angle  of  attack  in  all  cases.  In  the  limit  of 
vanishing  aspect  ratio,  the  present  theory  for  the  single  rectangular 
element  does  yield  values  of  9  which  appear  to  be  converging 
toward  Bollay's  value  of  half  the  angle  of  attack. 
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6.2  Triangular  Wing 

Calculations  have  been  carried  out  for  the  special  class  of 
swept  wings  of  triangular  planform  (i.e.,  delta  wings)  using  the 
method  of  Section  5.  Initially,  the  calculations  were  made  with 
a  double  iteration;  that  is,  with  iteration  on  both  the  loading 
and  the  shedding  angle  of  each  wing  element.  However,  for  reasons 
which  will  be  discussed  below,  the  final  calculations  were  carried 
out  by  iterating  on  the  loading  but  leaving  the  shedding  angle  6 
fixed.  Therefore,  convergence  of  the  swept  wing  solution  was 
investigated  both  with  and  without  iteration  on  8. 

6.2.1  Convergence  with  number  of  iterations 

Perhaps  the  key  to  the  success  or  failure  of  the  procedure 
described  in  Section  5.5  lies  in  the  rate  of  convergence  with 
iteration  on  the  loading  of  each  wing  element.  Therefore,  this 
was  investigated  first  using  four  rectangular  wing  elements. 
Furthermore,  since  the  rate  of  convergence  is  expected  to  deteri¬ 
orate  with  increasing  aspect  ratio,  the  highest  aspect  ratio 
(A  >  4)  to  be  employed  here  was  used  for  the  convergence  investi¬ 
gation.  Figure  8  shows  the  convergence  of  the  shed  vorticity  r*, 
the  first  loading  coefficient  7^*  (which  determines  C^  and 
x/c) ,  and  the  normal  force  coefficient  C„  with  number  of  itera- 
tions  on  the  loading.  It  can  be  seen  that  the  convergence  is 
generally  oscillatory  in  nature  and  that  in  this  case  six  iterations 
are  adequate  for  the  load  distribution,  whereas  only  three  itera¬ 
tions  are  required  to  converge  on  the  normal  force. 

Convergence  with  8  iterations  (which  are  carried  out  after 
convergence  on  7^,*)  is  illustrated  in  Figure  9.  It  can  be  seen 
that  convergence  is  again  oscillatory  and  that  four  iterations 
on  8  suffice  here  to  determine  8  within  about  0.25°.  On  the 
other  hand,  about  six  iterations  are  required  for  convergence  on 
the  normal  force. 
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6.2.2  Shedding  angle  calculations 

A  significant  difference  between  the  rectangular  and  tri¬ 
angular  wing  calculations  lies  in  the  fact  that  the  present  theory 
predicts  a  single  shedding  angle  for  the  rectangular  wing,  as  well 
as  a  single  value  for  the  shed  vorticity,  which  is  assumed  constant 
over  the  chord,  if  a  single  wing  element  is  used.  For  the  delta 
wing,  however,  the  method  yields  a  chordwise  distribution  of  these 
quantities  by  the  use  of  many  rectangular  wing  elements.  The  cal¬ 
culated  distribution  of  shedding  angle  6  over  the  chord  of  a 
delta  wing  of  aspect  ratio  4  at  20°  angle  of  attack  is  shown  in 
Figure  10  using  from  two  to  eight  wing  elements  to  represent  the 
delta  wing.  Here  we  see  again,  as  in  the  case  of  the  rectangular 
wing,  that  the  predicted  values  of  6  are  unrealistically  high. 

(in  this  case,  they  are  all  larger  than  the  angle  of  attack.)  But 
even  more  disturbing  is  the  observation  that  9  is  increasing  as 
the  number  of  wing  elements  increases.  This  may  be  associated 
with  the  increase  in  the  predicted  value  of  6  with  aspect  ratio, 
since  the  aspect  ratio  of  the  trailing  element  increases  as  the 
nunber  of  elements  is  increased. 

6.2.3  Convergence  with  number  of  harmonics 

The  convergence  of  the  normal  force  and  center  of  pressure 
on  the  delta  wing  was  next  investigated  using  various  numbers  of 
rectangular  elements  to  represent  the  delta  wings.  These  results 
are  shown  in  Figure  11,  and  it  can  be  seen  that  convergence  is 
quite  satisfactory  if  only  two  wing  elements  are  employed.  However, 
if  larger  nunbers  of  wing  elements  are  employed,  not  only  are  more 
harmonics  required  (as  would  be  expected  because  of  the  increased 
aspect  ratio  of  the  elements),  but  also  rather  large  bumps  are 
noted  in  the  variation  with  harmonics.  This  is  believed  to  be 
caused  by  the  fact  that  there  is  a  finite  discontinuity  in  the 
downwash  at  the  side  edge  of  each  rectangular  element.  Thus,  with 
certain  numbers  of  harmonics,  certain  spanwise  points  used  in  the 
Fourier  analysis  fall  close  to  these  discontinuities  and  cause 
apparent  convergence  difficulties.  In  order  to  check  this  hypothesis, 
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detailed  downwash  calculations  were  carried  out  for  one  case,  and 
the  fit  obtained  with  the  Fourier  series  was  investigated.  It  was 
found  that  by  adjusting  the  downwash  at  the  point  nearest  the  dis¬ 
continuity,  the  bumps  in  the  curves  of  Figure  11  could  be  reduced. 
It  is  also  felt  that  with  still  larger  numbers  of  wing  elements 
the  magnitude  of  each  bump  may  be  reduced  because  of  the  decrease 
in  the  magnitude  of  the  discontinuity. 

Perhaps  the  most  significant  results  of  Figure  11  are  that 

(a)  the  rate  of  convergence  with  number  of  wing  elements  is  slow, 

(b)  the  values  of  are  unrealistically  high,  and  (c)  the  value 

of  is  increasing  with  nvunber  of  wing  elements.  Inasmuch  as 

these  difficulties  with  C„  are  believed  to  be  associated  with 

N 

the  unrealistically  high  values  of  9  and  their  increase  with 
increasing  number  of  elements  (see  Fig.  10) ,  it  was  decided  to 
investigate  the  behavior  of  the  delta  wing  solution  without 
iterating  on  9 .  The  results  of  these  calculations  are  shown  in 
Figure  12,  and  it  can  be  seen  that  (a)  the  convergence  of  and 

x/c^  with  number  of  wing  elements  is  improved,  and  (b)  the  value 
of  Cjj  is  reduced,  even  for  the  rather  large  value  of  1.0  assumed 
for  the  ratio  of  0/a.  It  is  concluded  that,  in  applying  the 
present  theory  to  swept  wings ,  one  should  leave  0/a  as  a  param¬ 
eter,  because  of  the  unrealistic  predictions  of  0  for  the  assumed 
model . 
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6.2.4  Detailed  loading  coefficients  for  a  specific  case 

In  order  to  give  a  better  insight  into  the  behavior  discussed 
above,  detailed  loading  coefficients  for  each  rectangular  wing 
element  are  presented  in  Table  IV  for  a  delta  wing  of  aspect  ratio 
2  “0*5)  at  20°  angle  of  attack  with  the  value  of  0/a 

fixed  at  0.5.  For  this  calculation,  8  wings  and  19  nonzero  har¬ 
monics  are  employed,  and  the  characteristics  of  each  wing  element 
are  listed  in  the  following  table. 


Wing 

element 

No. 

!i 

s 

o 

°i 

°o 

C.L. 

°o 

L.L. 

°o 

0^ 

(deg) 

1 

0.0625 

0.5 

0.09375 

0.03125 

10 

2 

.1875 

1.5 

.21875 

.15625 

10 

3 

.3125 

.125 

2.5 

.34375 

.28125 

10 

4 

.4375 

.125 

3.5 

.46875 

.40625 

10 

5 

.5625 

.125 

4.5 

.59375 

.53125 

10 

6 

.6875 

.125 

5.5 

.71875 

.65625 

10 

7 

.8125 

.125 

6.5 

.84375 

.78125 

10 

8 

.9375 

,125 

7.5 

.96875 

.90625 

10 

where  C.L.  is  the  distance  from  the  apex  of  the  delta  wing  to 
the  control  line  of  each  element,  and  L.L.  is  the  distance  from 
the  apex  to  the  corresponding  lifting  line. 

The  calculated  loading  coefficients  on  each  element  are  listed 

in  Table  IV,  and  it  will  be  recalled  that  the  fourth  element  (Table 

IV(d))  has  the  same  geometry  as  the  isolated  rectangular  wing 

whose  loading  coefficients  are  listed  in  Table  III.  But,  before 

comparing  the  results  of  Tables  III  and  IV (d),  it  is  well  to  refer 

back  to  Equations  (75)  representing  the  delta  wing  problem.  It 

will  be  noted  that  all  of  the  mutual  interference  eunong  the  wing 

elements  is  contained  in  the  (b^),  and,  therefore,  appears  only 

J  ik 
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in  the  unseparated  solution 

r^*  ■  0.  Hence,  the  separated  solution  ^Yng* j  unaffected  by 

interference  and  should  be  the  same  as  for  an  isolated  rectangular 
wing  element  of  the  same  aspect  ratio.  Comparison  of  Tables  III 
and  IV(d)  shows  that  this  is  indeed  the  case.^° 

A  comparison  of  the  values  of  from  the  same  tables 

shows  that  in  the  case  of  the  delta  wing  (Table  IV (d))  the  mutual 
interference  among  the  wing  elements  severely  curtails  the  rate  of 
convergence;  that  is,  the  highest  harmonics  are  nearly  the  same 
order  of  magnitude  as  the  first  harmonic.  In  the  case  of  the 
corresponding  isolated  wing  element  (Table  III) ,  these  differed 
by  six  or  seven  orders  of  magnitude.  Looking  at  the  other  wing 
elements  of  the  delta  wing  (Table  IV) ,  we  see  that  this  difficulty 
is  typical  of  all  the  wing  elements  except  the  leading  element, 
vdiich  shows  convergence  similar  to  that  for  an  isolated  rectangular 
wing.  The  apparent  reason  for  this  difference  lies  in  the  fact 
that  the  present  mathematical  model  produces  finite  discontinuities 
in  the  downwash  for  all  wing  elements  except  the  first. 

Despite  the  above  difficulty,  which  may  be  expected  to  mani¬ 
fest  itself  in  the  calculation  of  such  detailed  quantities  as  span 
loading  and  downwash,  it  has  been  shown  in  the  previous  section 
and  in  Figure  12  that  the  gross  quantities  (C^^  and  for 

the  delta  wing  have  nevertheless  converged  with  8  wings  and  19  har¬ 
monics.  For  the  calculation  of  these  quantities  and  of  the  span 
loading,  we  require  the  values  of  shed  vorticity  P.*  for  each 
wing  element.  These  are  calculated  from  Equation  (51)  and  are 
listed  below. 


G'-a*), 


which 


is  obtained  by  setting 


^‘^he  small  differences  are  attributed  to  the  fact  that  the  numbers 
in  Table  III  were  calculated  on  an  IBM  1620  whereas  those  of 
Table  IV  were  calculated  on  an  IBM  7094. 
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Wing 

element 

^i* 

Mo. 

1 

0.34311 

2 

.17177 

3 

.11823 

4 

.09146 

5 

.07197 

6 

.05901 

7 

.04707 

8 

.03732 

The  total  normal  force  coefficient  and  center  of  pressure  for 
the  delta  wing  are  given  by  Equations  (87)  and  (90) .  The  result¬ 
ing  values  for  the  present  case  are 

C„  -  0.9292 

iM 

-  0.6272 

The  span  loading  on  the  delta  wing  can  be  calculated  from 

Equation  (94) .  The  resulting  span  load  distribution  is  tabulated 

as  a  function  of  y/s  in  Table  V.  It  will  be  noted  that  the 

o 

span  loading  goes  to  zero  at  ■  0.9375  which  corresponds  to 

the  side  edge  of  the  eighth  wing  element. 


-77- 


6.2.5  Comparisons  with  experiment 

The  calculated  variations  of  normal  force  and  center  of  pres¬ 
sure  with  angle  of  attack  for  delta  wings  of  aspect  ratio  1,  2,  3, 
and  4  are  compared  with  experiment  and  with  linear  theory  in 
Figures  13  through  16.  As  explained  earlier,  it  was  found  in  the 
course  of  the  calculations  that  iteration  on  the  shedding  angle  of 
the  separated  vortices  produced  unrealistically  high  values  of 
both  shedding  angle  and  normal  force  for  delta  wings.  Therefore, 
for  the  calculations  presented  in  Figures  13  through  16,  the  shed¬ 
ding  angle  was  treated  as  a  parameter  and  held  fixed  in  each  cal¬ 
culation. 

It  can  be  seen  from  Figures  13  through  16  that,  for  each 
aspect  ratio  investigated,  there  is  a  single  value  of  the  shedding 
angle  ratio  6/a  which  yields  good  agreement  with  the  experimental 
normal  force  for  all  angles  of  attack  up  to  the  experimental  fall¬ 
ing  off  of  the  lift-curve  slope.  This  latter  phenomenon  is  not 
predicted  by  the  present  theory  and  is  evidently  associated  with 
a  change  in  the  type  of  flow  separation  (as  in  the  case  of  rectan¬ 
gular  wings) .  The  calculated  normal  force  is  seen  to  increase 
with  e/a  in  all  cases.  (This  effect  is  shown  in  Figure  17  for 
various  aspect  ratios.)  On  the  other  hand,  the  center  of  pressure 
is  evidently  insensitive  to  d/a  and  agrees  well  with  experiment 
for  the  cases  where  data  are  available  (A  <■  1,  2,  4),  particularly 
at  the  higher  angles  of  attack.  It  is  noted  that  linear  theory 
also  shows  good  agreement  with  experiment  with  regard  to  center- 
of-pressure  location. 

It  is  evident  from  Figures  13  through  16  that  the  present 
theory  in  all  cases  yields  results  which  are  considerably  better 
than  linear  theory,  provided  that  the  appropriate  value  of  the 
ratio  d/a  is  selected  and  held  fixed.  The  required  6/a  is  a 
function  of  aspect  ratio,  and  the  resulting  variation  with  aspect 
ratio  is  presented  in  Figure  18. 
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Since  some  limited  experimental  values  of  the  shedding  angle 
0^  of  the  vortex  cores  over  slender  delta  wings  are  available  for 
aspect  ratios  of  1  and  1.67  (Refs.  24  and  25),  these  are  also 
shown  in  Figure  18.  One  is  again  cautioned  that  these  measure¬ 
ments  do  not  actually  correspond  to  experimental  values  of  0, 
since  the  separated  vortices  of  the  present  mathematical  model 
are  not  permitted  to  roll  up  into  a  single  core.  In  keeping  with 
this  fact,  it  is  seen  from  Figure  18  that  the  rolled-up  core  forms 
a  much  smaller  angle  with  the  plane  of  the  wing  than  do  the  separ¬ 
ated  vortices  of  the  assumed  mathematical  model  which  is  required 
to  give  the  experimental  normal  force. 

It  should  be  pointed  out  that  all  of  the  above  calculations 
were  performed  using  8  wing  elements  and  19  harmonics.  For  the 
highest  aspect  ratio  (4.0),  which  corresponds  to  the  lowest  shed¬ 
ding  angle  (0/a  *  0.1),  it  is  not  certain  that  the  mathematical 
solution  has  converged,  particularly  at  the  lowest  angle  of  attack. 

As  a  point  of  interest,  one  calculation  of  span  load  distri¬ 
bution  was  carried  out  for  an  aspect  ratio  2  delta  wing  at  20° 
angle  of  attack,  since  data  are  available  in  Reference  26.  The 
results  are  shown  in  Figure  19,  and  it  can  be  seen  that  discon¬ 
tinuities  are  indicated  in  the  theoretical  span  loading  at  the 
side  edge  of  each  rectangular  wing  element  employed  to  represent 
the  delta  wing.  It  is  observed  that  approximate  agreement  with 
experiment  would  be  obtained  by  fairing  out  the  discontinuities, 
but  it  is  evident  that  a  larger  number  of  wing  elements  (and 
harmonics)  would  be  required  for  accurate  prediction  of  the  span 
loading  on  delta  wings. 
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7.  CONCLUSIONS 

A  theoretical  analysis  has  been  developed  for  calculating  the 
low-speed  aerodynamic  loads  on  rectangular  wings  with  side-edge 
separation.  This  theory,  which  includes  lifting-line  theory  for 
the  case  of  no  separation,  yields  the  downwash  everywhere  in  the 
plane  of  the  wing  and  also  incorporates  the  Kutta  condition  of 
finite  velocity  along  the  side  edges.  Because  of  this  feature, 
it  has  been  possible  to  extend  the  theory  to  other  planforms  with 
sweptback  leading  edges.  An  iterative  technique  has  therefore 
been  devised  for  calculating  the  aerodynamic  characteristics  of 
sweptback  wings  with  leading-edge  separation  by  representing  the 
wings  as  a  system  of  elementary  rectangular  wings. 

Calculations  have  been  carried  out  for  both  rectangular  and 
delta  wings  of  various  aspect  ratios,  and  the  calculated  normal 
force,  center  of  pressure,  load  distribution,  and  shedding  angles 
compared  with  experimental  data.  The  following  conclusions  are 
drawn  from  the  calculated  results. 

7.1  Rectangular  Wings 

(1)  The  method  appears  to  converge  rapidly  with  respect  to 
number  of  harmonics,  except  at  low  angles  of  attack  and/or  high 
aspect  ratios.  These  cases,  however,  are  well  covered  by  lifting¬ 
line  theory  and  slender  wing  theory. 

(2)  A  characteristic  feature  of  the  spanwise  load  distribu¬ 
tion  for  a  rectangular  wing  with  side-edge  separation  is  found  to 
be  a  zero  slope  at  the  side  edges. 

(3)  For  aspect  ratios  from  zero  to  five,  the  normal  force 
predicted  by  the  present  theory  using  a  single  wing  element  agrees 
with  experiment  within  10  percent  up  to  the  onset  of  leading  edge 
stall,  which  is  not  included  in  the  theoretical  model. 

(4)  The  present  theory  evidently  predicts  unrealistically 
high  values  of  the  shedding  angle  for  the  separated  vortices, 
particularly  for  high  aspect  ratios  and  low  angles  of  attack. 
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(5)  The  predicted  normal  force  increases  with  shedding  angle 
of  the  separated  vortices  but  becomes  less  sensitive  to  shedding 
angle  as  the  aspect  ratio  increases. 

7.2  Delta  Wings 

(1)  For  delta  wings,  the  unrealistically  high  shedding 
angles  evidently  produce  correspondingly  high  predicted  values 
of  normal  force. 

(2)  If  the  shedding  angle  is  held  fixed,  the  present  theory 
is  capable  of  giving  substantial  improvement  over  other  theories 
in  predicting  both  normal  force  and  center  of  pressure  for  aspect 
ratios  from  one  to  four  over  a  wide  range  of  angle  of  attack, 
provided  an  appropriate  shedding  angle  is  selected.  A  plot  of 
the  required  shedding  angle  has  been  developed. 

(3)  The  predicted  center  of  pressure  is  relatively  insensi¬ 
tive  to  the  shedding  angle  of  the  separated  vortices. 

(4)  The  method  appears  to  converge  with  respect  to  number 
of  wing  elements  and  harmonics,  whether  or  not  one  solves  for  the 
shedding  angle  by  iteration.  For  aspect  ratios  up  to  four,  it 
appears  that  8  wings  and  19  nonzero  harmonics  suffice  for  2  percent 
precision  on  normal  force,  provided  that  the  separated  vortices 

lie  at  least  3^  above  the  wing  surface. 

(5)  Convergence  of  the  method  deteriorates  as  the  shedding 
angle  (and  consequently  the  angle  of  attack)  decreases  and  as  the 
aspect  ratio  increases .  That  is ,  larger  numbers  of  harmonics  are 
required  for  convergence  at  low  angles  of  attack  and/or  high 
aspect  ratios.  However,  linear  theory  suffices  for  those  cases. 

(6)  The  present  theory  affords  prediction  of  the  detailed 
span  loading  as  well  as  the  downwash  everywhere  in  the  flow  field. 
However,  accurate  prediction  of  the  span  loading  (and,  hence,  of 
the  downwash)  evidently  requires  larger  numbers  of  wing  elements 
and  harmonics  than  does  accurate  prediction  of  the  normal  force 
and  center  of  pressure. 


( 
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(7)  The  span  loading  on  a  swept  wing,  as  predicted  by  the 
present  method  using  8  wings  and  19  harmonics ,  exhibits  sizeable 
fictitious  discontinuities  at  the  side  edge  of  each  wing  element. 
However,  when  these  were  faired  out,  the  resulting  span  load  dis¬ 
tribution  compared  well  with  experiment  for  the  one  case  which 
was  calculated  for  an  aspect  ratio  two  delta  wing  at  20°  angle  of 
attack. 

(8)  The  present  theory  evidently  does  not  predict  the  fall¬ 
ing  off  of  lift-curve  slope  near  the  stall.  Thus,  for  delta  wings 
of  aspect  ratio  three  and  higher,  the  theoretical  normal  force  is 
too  high  for  angles  of  attack  above  20°. 
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APPENDIX  A 


UPWASH  INDUCED  BY  A  HORSESHOE  VORTEX 
LYING  ABOVE  THE  WING 


To  determine  the  upwash  w  induced  at  an  arbitrary  point 
(x,y,z)  by  a  horseshoe  vortex  which  is  inclined  to  the  wing 
surface,  let  us  consider  the  following  sketch: 


z 


Here  we  have  a  rectangular  wing  of  semispan  s  lying  in  the  x,y 
plane  and  a  horseshoe  vortex  of  strength  r.  This  vortex  consists 
of  three  parts,  a  bound  part,  of  semispan  s,  lying  in  the  plane 
of  the  wing  at  a  distance  ^  behind  the  leading  edge,  and  a  right 
and  left  trailing  vortex  each  of  which  is  inclined  at  an  angle  d 
to  the  wing  as  shown.  Let  us  divide  the  upwash  w  produced  at 
point  (x,y,z)  into  the  components  produced  by  each  part  of  the 
horseshoe  vortex.  That  is,  let 

w(x,y,z)  -  Wg(x,y,z)  +  Wj^(x,y,z)  +  Wj^(x,y,z)  (A.l) 


The  upwash  induced  by  the  bound  vortex  can  be  determined 
with  the  aid  of  the  following  sketch: 


Thus,  the  total  velocity  induced  at  point  P  by  the  bound  vortex 
of  span  2s  is  (Ref.  27) 

“  4^7  “2^  \ 

where  ia  is  the  unit  vector  in  the  direction  of  q_.  Hence 
the  upwash  in  the  positive  z  direction  is  given  by  (see  sketch) 


cos  <p 


r..gos,  Sl> 

4Trr 


(cos  a  +  cos  a  ) 
1  2 


(A. 3) 


where  T  is  the  unit  vector  in  the  z  direction. 


The  preceding  sketch  yields  the  following  relationships 


- 


cos  a  — - - 

V (x  -  4) ^  +  (s  +  y) ® 

cos  ■  — - -  y.. - - 

V  (x  -  i)^  +  (s  -  y)2 

so  that  substitution  into  Equation  (A. 3)  gives 


-  _ 1C.(2C  _ 

s  +  y 

®  47r 

^TT 

^(x  -  +  z  J 

■V  (x  -  +  (s  +  y)  +  z* 

(x  -  +  (s  -  y)^  +  z* 


(A. 5) 


The  next  component  of  the  upwash  at  point  P  to  be  determined 
is  that  due  to  the  right  trailing  vortex.  For  this  purpose  let 
us  look  at  the  next  sketch. 

z  s  y 


A-4 


If  we  apply  the  relationship  given  by  Equation  (A. 2)  to  this  case, 
we  have 


>  (cos  a  +  cos  a  )  i 
4vr  '3  4  q» 


(A. 6) 


where  iq^  is  the  unit  vector.  Since  the  vortex  we  are  dealing 
with  is  semi-infinite,  approaches  zero  degrees  and 


cos  a  “1.0 
4 


(A. 7) 


Let  us  now  write  the  expressions  for  two  unit  vectors,  one  in  the 
direction  R  to  Q  and  the  other  in  the  direction  R  to  P.  These 
are 

T„_  «  cos  6  T  +  sin  0  T  (A. 8) 
£\Vt  X  z 


and 


(x  -  -  (s  -  y)Ty  +  zT^ 

V  (x  -  +  (s  -  y)*  +  z^ 


(A. 9) 


From  these  two  unit  vectors  we  find 

.  _ (x  -  ^)cos  e  +  z  sin  0 


cos 


V 


(A.  10) 


(x  -  +  (s  -  y)®  +  z® 


r  ■  |^p|  sin  a  ■  j  ^P  1  -  cos® 


and 


^(s  -  y)^  +  ^(x  -  |)sin  0  -  z  cos  0J  ^ 


.  ^RQ  ^  '•RP  . 

‘Jr  ‘*3  ^ 


1/2 


(A. 11) 
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This  cross  product  is,  from  Equations  (A. 8)  and  (A. 9) 

IT  X  ftp  c 

■  - “  r  1  ~  y)8in  6  +  ^(x  -  ^)sin  6  -  z  cos  sj  T 


-  (s  -  y)cos  6  T, 


y 

(A. 12) 


Therefore,  substitution  of  Equations  (A. 7),  (A. 10),  (A. 11),  and 

(A. 12)  into  Equation  (A. 6)  will  give  the  expression  for  'q„.  Thus, 

for  the  upwash  induced  at  P  (i.e.,  the  component  of  q^  in  the 

T  direction) ,  we  find 
z  ’ 

%  ‘  ^z 


.  -f 

4Tr 


is,  z..i)s9B_e_ 


(s  -  y)^  +  [^(x  -  4) sin  0  -  z  cos  ej’ 


(x  -  £)co8  6  -f  z  sin  6  ^ 

V(x  -  +  (s  -  y)^  +  z* 


(A. 13) 


From  symmetry  considerations,  the  ejqpression  for  the  upwash 
induced  at  point  P  by  the  left  trailing  vortex  can  immediately 
be  written  as 


w. 


-C-  J 

r 

[s  +  y)cos  0 

4ir  1 

1  (s  +  y)®  + 

(x  -  ^)8in  0  -  z  cos  0^ 

V 


(x  -  ^)coB  6  +  z  sin  0 
(x  -  ^)  +  (s  +  y)  +  z 


+  1 


(A.  14) 


Finally,  the  total  upwash  w  due  to  the  horseshoe  vortex 
of  strength  F  shown  in  the  first  sketch  of  this  appendix  is 


W  -  Wb  +  W^  +  Wj^ 


where 


r(x  -  e) 


r 


a  +  V 


APPENDIX  B 


INFLUENCE  OF  THE  GROUND  PLANE 

If  the  wing  in  question  is  in  the  proximity  of  a  ground 
plane,  then  the  upwash  at  the  wing  surface  will  be  altered  in  the 
same  fashion  as  if  the  ground  plane  were  replaced  by  an  inverted 
"image"  wing  located  as  shown  in  the  sketch  below  (Ref.  28): 


r 


B-2 


Thus,  we  can  rewrite  the  boundary  condition  of  no  flow  through  the 
wing  surface  (Eq.  (9) )  in  the  form 

(B.l) 

where  the  subscript  g  denotes  the  contribution  of  the  ground 
plane  or  image  system. 

It  can  be  seen  from  the  above  sketch  that  the  z  component 
of  the  velocity  produced  at  point  P(x,  y)  on  the  wing  by  the 
image  system  is  equal  to  the  z*  component  of  velocity  produced 
at  the  image  point  P'(x',  y')  by  the  actual  vortex  system  asso¬ 
ciated  with  the  real  wing.  Therefore,  we  can  determine  the  upwash 

w  *(x,  y)  and  w  (x,  y)  of  Equation  (B.l)  by  writing  expressions 
y  y 

for  the  velocity  component  in  the  z'  direction  produced  at  the 
image  point  P'  by  the  actual  wing.  For  this  purpose,  it  can  be 
seen  from  the  sketch  below  that  we  shall  require  both  the  x  and 
z  components  of  velocity  induced  at  the  point  P’. 


Ground  plane 


That  is , 


w. 


Mi  •  m  -V  cos  2a  +  u  sin  2a 
9 


(B.2) 
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where  w  and  u  are  the  z  and  x  components  of  velocity 
induced  at  P' . 

The  coordinates  x^,  z^  (in  the  x,  y,  z  system)  of 

the  image  point  P'  are  required  for  use  in  the  general  expressions 
for  w  and  u.  But,  since  the  y  coordinate  of  the  image  point 
is  the  same  as  that  for  the  corresponding  point  in  the  wing  itself, 
we  have 

y^  -  y  ^  (B.3) 

The  coordinates  x^^  and  z^  can  be  obtained  from  the  following 
sketch,  which  is  a  cut  through  the  y  “  y^  plane. 


-i-  (c  -  x)sin  a  cos  a  (B.4) 


+  2  [h,pg  +  (c  -  x)sin  a]  sin  a 


Thus,  we  see  that 

*i  “  ■  ^  [Nt: 

and 

x^  «  X  -  Zji^  tan  a  ■  X 


(B.5) 
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Having  the  coordinates  x^,  of  the  image  point  P', 

as  given  by  Equations  (B.3),  (B.4),  and  (B.5),  we  can  now  express 
the  u  and  w  velocity  components  at  that  poxnt  by  using  the 
results  of  Appendix  A.  Thus,  the  w  component  is  given  directly 
by  Equations  (A. 15)  if  we  replace  x  by  x^  and  z  by  z^. 
Similarly,  the  u  components  at  point  P{x,  y,  z)  are  obtained 
from  Equations  (A. 2)  and  (A. 6)  by  writing  the  appropriate  dot 
products  of  the  velocity  vectors.  Thus,  from  Equation  (A. 2), 
we  have  (see  sketch,  p.  A-2) 


sin  A  ■  ^  ^  (cos  a  +  cos  a  ) 

4'irr  1  2 


(B.6) 


Therefore,  using  the  relationships  of  Equation  (A,4) ,  we  find  that 


_ _ LS _ 

r _ s  +  y _ 

B  . 

47r  1 

4.  - 

_(x  - 

_ S- -■  i 

r 

1 

^  +  (s  +  y)^  +  z^ 

V (x  -  ^)^  +  (s  -  y)^  +  z^_ 

(B.7) 


By  the  same  procedure,  we  find  from  Equations  (A. 6)  through 
(A. 12)  that 


.  JL 


c- 


(s  -  y)sin  e 


^  (s  -  y)^  +  [(x  -  4)sin  6  -  z  cos 
(x  -  ^)cos  9  +  z  sin  e 


) 


V (x  - 


+  1 


?)^  +  (s  -  y)^  +  z® 


(B.8) 


B-5 


From  symmetry  considerations,  the  corresponding  expression  for 
induced  at  P  by  the  left  trailing  vortex  can  immediately 
be  written  as 

u  .  _L  r _ (s  +  9 _ 

L  ^  (s  +  y)®  +  ^(x  -  ^)8in  6  -  z  cos 

(x  -  ^)cos  0  +  z  sin  0  ^ 

V(x  -  +  (s  +  y)^  +  2® 

Finally,  the  u  and  w  conqponents  induced  at  the  image 
point  P' (x^,  y^,  are  obtained  directly  from  Equations  (A. 15), 

(B.7),  (B.8),  and  (B.9)  by  simply  replacing  x  and  z  by  x^^ 
and  z^  as  given  by  Equations  (B.4)  and  (B.5).  Then  the  actual 
upwaah  Wg  produced  at  the  point  P  by  the  image  wing  (i.e.,  by 
the  ground  plane)  is  given  by  Equation  (B.2). 

It  must  be  recalled  that  the  above  expressions  give  only  the 
velocities  induced  by  a  single  vortex,  whereas  the  additional 
upwash  components  required  in  the  boundary  condition  of  Equation  (B.l] 
are  those  produced  by  the  entire  vortex  system  associated  with  the 
image  wing.  Therefore,  in  order  to  obtain  w  *  and  w  ,  we  must 
integrate  the  above  expressions  in  a  manner  analogous  to  that 
employed  in  the  main  body  of  the  present  report  (see  Sections  4.2.2 
and  4.2.3).  The  integrations  required  will,  of  course,  be  more 
difficult  than  those  of  the  main  analysis,  since  they  now  involve 

points  lying  out  of  the  plane  of  the  wing  (i.e.,  z  ^  0) . 

Inasmuch  as  the  detailed  analysis  and  programing  of  the 

influence  of  the  ground  plane  represents  a  formidable  task,  it  is 

well  to  re-excunine  the  validity  of  the  assumptions  made  in  this 
analysis  before  embarking  on  such  an  endeavor.  That  is,  one 
should  perhaps  first  consider  whether  the  mathematical  model  is 
sufficiently  realistic  to  justify  this  extension.  In  particular, 
it  is  pointed  out  that  the  shed  vortices  associated  with  the 
lifting  line  are  assumed  to  extend  in  the  chordwise  direction 


B-6 


while  the  separated  vortices  are  assumed  to  be  shed  at  an  angle  6 
above  the  wing.  Even  in  the  absence  of  a  ground  plane,  this  assump¬ 
tion  is  of  course  an  approximation  to  the  real  vortex  system  in 
which  the  vortices  roll  up  into  two  concentrated  cores.  But  as 
the  wing  approaches  the  ground,  an  additional  complication  arises 
which  is  not  accounted  for  in  the  present  analysis.  The  shed 
vortices  are  deflected  upward  by  the  influence  of  the  ground  plane 
and  are  also  driven  outward  in  the  manner  of  a  vortex  ring  which 
expands  when  approaching  a  flat  surface.  This  latter  effect  has 
been  observed  experimentally  behind  a  delta  wing  model .  These 
effects  will,  of  course,  become  increasingly  important  as  the 
ground  plane  is  approached,  but  their  quantitative  effect  upon 
the  calculated  results  is  not  known. 


TABLE  I.-  DOWNWASH  IN  WAKE  BEHIND  ELLIPTICALLY  LOADED 

LIFTING  LINE. 
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TABLE  II.-  VARIATION  IN  RECTANGULAR  WING  CHARACTERISTICS 
WITH  NUMBER  OF  NON-ZERO  FOURIER  HARMONICS.  (CONTROL 
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TABLE  II.-  CONCLUDED. 
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TABLE  III.-  CALCULATED  LOADING  COEFFICIENTS  FOR  AN 
ISOLATED  RECTANGULAR  WING  ELEMENT  OF  ASPECT  RATIO  3.5 
(a  -  20°,  e  -  10°,  -  19,  x^c  -  0.75). 
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7.7604X10"“ 

4.8911X10"® 

23 

1.8738X10"® 

2.9531X10"“ 

1.8645X10"® 

25 

2.5503X10"® 

5.3567X10"® 

3.3813X10"® 

27 

-1.1670X10"® 

-3.4224X10"® 

-2.1704X10"® 

29 

-6.4205x10"® 

-4.7103X10"® 

-2.9774x10"® 

31 

-3.3350X10"® 

-3.4346X10"® 

-2.1697X10"® 

33 

7.2535X10"® 

-1.8373x10"® 

-1.1516X10"® 

35 

-1.4820x10"® 

-2.3315x10"® 

-1.6188X10"^ 

37 

-7.7647X10"^^ 

2.8269x10"® 

1.7823X10"^ 

TABLE  IV.-  CALCULATED  LOADING  COEFFICIENTS  FOR  EACH 
RECTANGULAR  WING  ELEMENT  OF  A  DELTA  WING  OF  ASPECT 
RATIO  2.0  (a  -  20®,  0^  -  10°,  (Xq/c)^  -  0.75,  -  19) . 


(a)  Wing  element  number  1. 


n 

Ms 

y  * 

1 

0.2155 

-0.5293 

0.3392X10"^ 

3 

-.1248X10*® 

-.9414x10*^ 

-.3242X10*^ 

5 

-.1106X10"* 

-.5724x10*® 

-.1975X10"® 

7 

.2877X10*® 

.1199X10*® 

.4117X10"® 

9 

-.2560X10*® 

.2273X10*® 

.7798X10"“ 

11 

-.2697X10*® 

-.1943X10"“ 

-.6668X10"® 

13 

.3009x10"® 

-.9041X10"® 

-.3099X10"® 

15 

.1324x10*® 

.  1016x1  O*"^ 

.3621x10"® 

17 

-.1151X10*® 

.3423x10*® 

.1163X10"® 

19 

.2171X10*® 

.1988X10*’^ 

.8995x10"® 

21 

-.6955X10*® 

-.2789x10*'^ 

-.1026X10"'^ 

23 

.3579X10*® 

.7980X10*® 

.3096X10"® 

25 

-.3787x10*® 

.9309x10*® 

.2815x10*® 

27 

.1605x10*® 

-.9566X10*® 

-.1676X10"® 

29 

.2140X10*® 

-.8673x10*® 

-.8358X10"® 

31 

-.1896X10*® 

.6237X10*® 

.2434X10"® 

33 

-.1955X10*® 

.8279x10*® 

.2645X10"® 

35 

-.8298x10*® 

.1249x10*® 

.3458X10"^ 

37 

.5229x10*® 

-.1171X10*® 

-.3496X10"’’ 

TABLE  IV.-  CONTINUED 


(b)  Wing  element  number  2 . 


n 

^n*)„ 

Y  * 

'n 

1 

0.1584 

-0.6940 

0.3924X10“^ 

3 

.6306X10"^ 

-.3015 

.1126X10“^ 

5 

-.4885x10"^ 

-.1012 

-.6625X10“^ 

7 

.3059X10"^ 

-.2411x10“^ 

.2645X10“^ 

9 

-.1599X10“^ 

-.1648X10“® 

-.1628X10“^ 

11 

.6325X10"^ 

.1864x10“® 

.6645X10“® 

13 

-.4304X10“® 

.1119x10“® 

-.2380X10“® 

15 

-.2920X10“® 

.3221X10“® 

-.2865X10“® 

17 

.4364X10“® 

.1190X10““ 

.4366X10“® 

1? 

-.4204x10“® 

-.3980X10““ 

-.4210X10“® 

21 

.2817X10“® 

-.2279x10““ 

.2814X10“® 

23 

-.8444X10“® 

-.6200x10“® 

-.8454x10“® 

25 

-.9381X10“® 

.1882x10“® 

-.9381X10“® 

27 

.1890X10“® 

.1051X10“® 

.1890x10“® 

29 

-.1766X10“® 

.5290x10“® 

-.1766X10“® 

31 

.7988X10“® 

.1458X10“® 

.7988X10“® 

33 

.4269X10“® 

.1400x10“'^ 

.4269X10“® 

35 

-.8806x10“® 

.2557X10“® 

-.8806X10“® 

37 

.3724X10“® 

-.3111x10“® 

.3723X10“® 

I 

i  TABLE  IV.-  CONTINUED. 


(c)  Wing  element  number  3. 


n 

Mu 

Ms 

y  * 

'n 

1 

0.9065x10"^ 

-0.7347 

0.3776x10"® 

3 

.1178 

-.4157 

.6868X10"^ 

5 

-.2504X10"^ 

-.1917 

-.4771X10"^ 

7 

-.4153X10"^ 

-.7680X10“^ 

-.5061X10"^ 

9 

.4409X10"^ 

-.2399X10“^ 

.4126X10“^ 

11 

-.1425X10"^ 

-.3903X10“® 

-.1471X10"^ 

13 

-.8608X10“® 

.1560X10"® 

-.8424X10"® 

15 

.1121X10"^ 

.1873X10“® 

.1143X10“^ 

17 

-.2247x10“® 

.1060X10“® 

-.2121X10“® 

19 

-.4651x10“® 

.3971X10“® 

-.4604X10"® 

21 

.3053x10“® 

.7052X10““* 

.3061X10“® 

23 

.3506X10“® 

-.3419x10"“* 

.3502x10"® 

25 

-.7061x10“® 

-.4214X10"“* 

-.7066x10“® 

27 

.3765x10“® 

-.2472X10"“* 

.3762X10“® 

29 

.2521X10“® 

-.9463X10“® 

.2520x10“® 

31 

-.5000X10“® 

-.1617X10"® 

-.5001x10“® 

33 

.1624X10“® 

.7652X10"® 

.1624x10"® 

35 

.1778X10“® 

.1922X10"® 

.1778X10“® 

37 

-.1139x10“® 

.3864X10"® 

-.1139x10"® 

TABLE  IV.-  CONTINUED 


(d)  Wing  element  nvutiber  4. 


n 

6'"*X 

y  * 

1 

0.6064X10"^ 

-0.7483 

-0.7796X10"® 

3 

.1080 

-.4894 

.6326X10"^ 

5 

.2722X10"^ 

-.2616 

.3293X10"® 

7 

5315X10"^ 

-.1285 

-.6490X10"^ 

9 

-.1071X10“^ 

-.5547X10"-^ 

-.1578X10"^ 

11 

.3937X10“^ 

-.1920X10"^ 

.3762X10"^ 

13 

-.4958X10"® 

-.3758X10"® 

-.5301X10"® 

15 

-.2755X10"^ 

.1322X10"® 

-.2743X10"^ 

17 

.1824x10"^ 

.2075X10"® 

.1843X10"^ 

19 

.5509X10"® 

.1486X10"® 

.5645x10"® 

21 

-.7850x10"® 

.7758X10"® 

-.7779X10"® 

23 

-.3288X10"® 

.2954X10"® 

-.3261X10"® 

25 

.1978x10"® 

.5360X10"“ 

.1983X10"® 

27 

.7397X10"® 

-.3428X10"“ 

.7394x10"® 

29 

-.3141X10"® 

-.4726X10"“ 

-.3146X10"® 

31 

-.9005x10"® 

-.3440x10"“ 

-.9008X10"® 

33 

.5349X10"® 

-.1831X10"“ 

.5347X10"® 

35 

.5580X10"® 

-.1887X10"® 

.5580X10"® 

37 

-.4149x10"® 

.2268X10"® 

-.4149X10"® 

TABLE  IV.-  COMTINUBD. 


(e)  Wing  element  number  5. 


n 

1 

0.3559x10"^ 

-0.7527 

-0.1859X10"^ 

3 

.8851X10”^ 

-.5412 

.4956X10"^ 

5 

.6267x10"^ 

-.3170 

.3985X10"^ 

7 

-.3505X10"-^ 

-.1743 

-.4760X10"^  1 

9 

-.4086X10'’^ 

-.8825X10"^ 

-.4721X10"^ 

11 

.1596X10"^ 

-.3933X10"^ 

.1312X10"^ 

13 

.3099X10"^ 

-.1406X10”^ 

.2998X10"^ 

15 

-.1223x10"^ 

-.2687X10"® 

-.1243X10"^ 

17 

-.2307x10“^ 

.1385X10"® 

-.2297X10"^ 

19 

.1712X10"^ 

.2132X10"® 

.1728X10"^ 

21 

.4498X10"® 

.1682X10"® 

.4619X10"® 

23 

-.6930X10"® 

.1015X10"® 

-.6857X10"® 

25 

.3880X10"® 

.4886x10"® 

.3915X10"® 

27 

-.8885X10"® 

.1688X10"® 

-.8873X10"® 

29 

.4527X10"® 

.1365X10"* 

.4528X10"® 

31 

.1052X10"”^ 

-.3891X10"* 

.1052X10"^ 

33 

-.9162X10"® 

-.3805X10"* 

-.9165X10"® 

35  1 

-.3712X10"® 

-.4074X10"* 

-.3715x10"® 

37 

.4023X10"® 

-.2812X10"* 

.4021X10"® 

1 

u 


TABLE  IV.-  CONTINUED. 


(f)  Wing  element  number  6. 


n 

O'"*! 

Yn* 

1 

0.2618X10“^ 

-0.7535 

-0.1827X10"^ 

3 

.7148X10"^ 

-.5794 

.3728X10"^ 

5 

.6104X10"^ 

-.3623 

.3966X10"^ 

7 

-.1285x10"® 

-.2143 

-.1393X10"^ 

9 

-.4052X10”^ 

-.1195 

-.4757X10"^ 

11 

-.2327X10"^ 

-.6118X10"^ 

-.2688X10"^ 

13 

.3458X10"^ 

-.2744X10"^ 

.3296X10"^ 

15 

.1702X10"^ 

-.9671X10"® 

.1645X10"^ 

17 

-.1744X10"^ 

-.1446X10"® 

-.1752X10"^ 

19 

-.1421X10"^ 

.1593X10"® 

-.1412X10"^ 

21 

.7034X10"® 

.2156x10"® 

.7161x10"® 

23 

. 8889X10"® 

.1756X10"® 

.8992X10"® 

25 

-.1244x10"® 

.1143x10"® 

-.1176X10"® 

27 

.2460X10"® 

.6236X10"® 

.2496X10"® 

29 

-.1128X10"^ 

.2782X10"® 

-.1127X10"^ 

31 

-.1129X10“® 

.9031X10"“ 

-.1124X10"® 

33 

.8922X10"® 

.1318X10"* 

.8923X10"® 

35 

.1978X10"® 

.1668X10"® 

.1988X10"® 

37 

-.4042X10"® 

.1667X10"® 

-.4032X10"® 

i 


I: 


TABLE  IV.-  CONTINUED. 


(g)  Wing  element  number  7. 


n 

(vX 

y  * 

'n 

1 

0.1504X10“^ 

-0.7527 

-0.2038X10"^ 

3 

.5024X10“^ 

-.6086 

.2159X10"^ 

5 

.6197X10“^ 

-.4002 

.4313X10"^ 

7 

.2096X10“^ 

-.2494 

.9227X10"® 

9 

-.3282X10"^ 

-.1487 

-.3982X10"^ 

11 

-.2544X10"^ 

-.8319X10"^ 

-.2936X10"^ 

13 

-.3194X10“^ 

-.4245X10"^ 

-.5192X10"® 

15 

.3308X10"^ 

-.1877X10"^ 

.3220X10"^ 

17 

.1888X10"® 

-.6201X10"® 

.1596X10"® 

19 

-.1562X10"^ 

-.3332X10"® 

-.1564X10"^ 

21 

-.4374X10"® 

.1826X10"® 

-.4288X10"® 

23 

.5223x10"® 

.2167X10"® 

.5325x10“® 

25 

.8419X10"® 

.1776X10"® 

.8503X10"® 

27 

-.1442X10"^ 

.1213X10"® 

-.1436X10"^ 

29 

.8713X10"® 

.7236X10"® 

.8747X10"® 

31 

.5125X10"® 

.3795X10"® 

.5143X10"® 

33 

-.6490X10"® 

.1690X10"® 

-.6482X10"® 

35 

-.3872X10"® 

-.1372X10"® 

-.3878X10"® 

37 

_ 

.3965X10"® 

-.1861X10"® 

.3956X10"® 

I 


I 

1  TABLE  IV.-  CONCLUDED. 


(h)  Wing  element  number  8. 


n 

(^"*)s 

y  * 

’n 

1 

0.7368X10"^ 

-0.7512 

-0.2067X10"^ 

3 

.4651X10"^ 

-.6314 

.2294x10"^ 

5 

.5005x10“^ 

-.4325 

.3390x10"^ 

7 

.1983x10"^ 

-.2804 

.9367X10"® 

9 

-.8729X10"^ 

-.1755 

-.1528X10"^ 

11 

-.3937x10"^ 

-.1046 

-.4328X10"^ 

13 

-.6698X10"® 

-.5815x10"^ 

-.8869X10"® 

15 

.1373x10"^ 

-.2925X10"^ 

.1264x10"^ 

17 

.3284X10"^ 

-.1247X10"^ 

.3237X10"^ 

19 

-.9760x10"® 

-.3561x10"® 

-.9892x10"® 

21 

-.2146X10"^ 

.5685X10"® 

-.2144x10"^ 

23 

-.5058x10"® 

.2032X10"® 

-.4982x10"® 

25 

.7683x10"® 

.2172X10"® 

.7764X10"® 

27 

.1418X10"^ 

.1772X10"® 

.1425x10"^ 

29 

-.4515x10"® 

.1246X10"® 

-.4469x10"® 

31 

.1036X10"® 

.7753X10"® 

.1065x10"® 

33 

-.1018x10"^ 

.4086x10"® 

-.1016x10"^ 

35 

.1435x10"® 

-.5035x10"* 

.1433x10"® 

37 

.3316x10"® 

-.1795X10"® 

.3310X10"® 

TABLE  V.-  CALCULATED  SPAN  LOAD  DISTRIBUTION  FOR 
A  DELTA  WING  OF  ASPECT  RATIO  2.0  AT  20°  ANGLE 
OP  ATTACK  (8  wings,  19  harmonics,  9/a  ■  0.5). 


y/Sq 

c,  c/c 

1  o 

y/s 

o 

=jO/C^ 

0.0000 

0.2072 

0.4625 

0.1122 

.0125 

.2073 

.4875 

.1159 

.0250 

.2074 

.5125 

.1190 

.0375 

.2075 

.5375 

.1224 

.0500 

.2077 

.5625 

.0854 

.0625 

.1864 

.5875 

.0898 

.0875 

.1868 

.6125 

.0934 

.1125 

.1874 

.6375 

.0972 

.1375 

.1880 

.6625 

.0993 

.1625 

.1893 

.6875 

.0608 

.1875 

.1578 

.7125 

.0650 

.2125 

.1594 

.7375 

.0673 

.2375 

.1624 

.7625 

.0711 

.2625 

.1653 

.7875 

.0763 

.2875 

.1675 

.8125 

.0421 

.3125 

.1330 

.8375 

.0448 

.3375 

.1366 

.8625 

.0439 

.3625 

.1406 

.8875 

.0420 

.3875 

.1431 

.9125 

.0383 

.4125 

.4375 

.1454 

.1080 

.9375 

.0000 

Aspect  ratio.  A 


Figure  2.-  Dependence  of  lift-curve  slope  on  location  of  control 
line  for  rectangular  wings,  from  Equation  (64). 
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(a)  Normal  force. 
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Figure  7.-  Concluded. 
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(b)  First  loading  coefficient. 


(c)  Delta  wing  normal  force  coefficient 
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Figure  10.-  Chordwise  variation  of  shedding  angle  on  aspect 
ratio  4  delta  wing  at  20*^  angle  of  attack. 


rectangular 

elements 


(a)  Delta  wing  normal  force  coefficient. 

Figure  11.-  Convergence  of  aerodynamic  characteristics  with  numbers 
of  rectangular  elements  and  nonzero  harmonics. 


(b)  Delta  wing  center  of  pressure 
■  Figure  11.—  Concluded. 
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Figure  12.-  Convergence  without  6  iteration 
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(a)  Normal  force. 

Figure  13.-  Aerodynamics  of  a  delta  wing  of  aspect  ratio  1  for 
various  shedding  angles  (8  wings,  19  harmonics). 


Symbol 


0.5  ^ 

.sj 


Experimental  data 
(Ref.  25) 


Figure  14.-  Aerodynamics  of  a  delta  wing  of  .aspect  ratio  2 
for  various  shedding  angles  (8  wings,  19  harmonics) . 


(a)  Normal  force. 

Figure  15.-  Aerodynamics  of  a  delta  wing  of  aspect  ratio  3 
for  various  shedding  angles  (8  wings,  19  harmonics). 
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Figure  15.-  Concluded. 


(a)  Normal  force. 


Figure  16.-  Aerodynamics  of  a  delta  wing  of  aspect  ratio  4 
for  various  shedding  angles  (8  wings,  19  harmonics). 
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Figure  18.-  Value  of  shedding  angle  required  for  accurate 
prediction  of  normal  force  on  delta  wings. 


APPROVED  DISTRIBUTION  LIST 


I 


I 

CONTRACT  Nonr 


'  Chief,  Bureau  of  Naval  Weapons 
{RAAD-3) 

I  Department  of  the  Navy 
I  Washington  25,  D.  C. 

j  Chief,  Bureau  of  Naval  Weapons 
1  (RAAD-33) 

Department  of  the  Navy 
j  Washington  25,  D.  C. 

Chief,  Bureau  of  Naval  Weapons 
(RAAD-34) 

Department  of  the  Navy 
Washington  25,  D.  C. 

I'  Chief,  Bureau  of  Naval  Weapons 
.  (RA-4) 

Department  of  the  Navy 
j- Washington  25,  D.  C. 

Chief,  Bureau  of  Naval  Weapons 
(R-55) 

Department  of  the  Navy 
Washington  25,  D.  C. 

I  Chief,  Bureau  of  Naval  Weapons 
I  (RR-25) 

Department  of  the  Navy 
I  Washington  25,  D.  C. 

Chief,  Bureau  of  Naval  Weapons 
(RRRE-4) 

I  Department  of  the  Navy 
'  Washington  25,  D.  C. 

1  Commanding  Officer  and  Director 
I  David  Taylor  Model  Basin 
Aerodynamics  Laboratory 
j  Washington  25,  D.  C. 

Chief  of  Naval  Research  (Code  461) 

.  Department  of  the  Navy 
I  Washington  25,  D.  C.  (6  copies) 

I 

I 


3103 (00) 


Chief  of  Naval  Research  (Code  438) 
Department  of  the  Navy 
Washington  25,  D.  C. 

Commanding  Officer 

Office  of  Naval  Research  Branch  Office 
Navy  #100,  Box  39,  F.  P.  O. 

New  York ,  New  York 

ATTN:  Head,  Documents  Section 

(2  copies) 

Commanding  Officer 

Office  of  Naval  Research  Branch  Office 

346  Broadway 

New  York  1 3 ,  New  York 

Commanding  Officer 
Office  of  Naval  Research  Branch 
86  E.  Randolph  Street 
Chicago  1,  Illinois 

Commanding  Officer 

Office  of  Naval  Research  Branch  Office 
1030  E.  Green  Street 
Padadena,  California 

Director 

Naval  Research  Laboratory 
Technical  Information  Office 
Washington  25,  D.  C.  (6  copies) 

Commander 

Army  Material  Command 
Department  of  the  Army 
Washington  25,  D.  C. 

ATTN:  AMCRD-RS-PE-A  (2  copies) 

Commanding  Officer 
U.  S.  Army  Transportation  Research 
Command 

Fort  Eustis,  Virginia 

ATTN:  SMOFE-TD  (1  copy) 

ATTN:  Research  Reference  Center 

(1  copy) 


r 


I  Army  Research  Center 

Physical  Sciences  Division 
3045  Columbia  Pike 
Arlington  4,  Virginia 
I  ATTN:  Mr.  R.  Ballard 

U.  S.  Air  Force  (SRGL) 
j  Office  of  Scientific  Research 
!  Washington  25,  D.  C. 

j  Aeronautical  Systems  Division  ■ 
j  Deputy  for  Technology 

Wright-Patterson  AFB,  Ohio 
attn:  asrsms 

I 

•  Aeronautical  Systems  Division 
Deputy  for  Technology 
Wright-Patterson  AFB,  Ohio 

_  attn:  asrssc 

I'  Aeronautical  Systems  Division 
Deputy  for  Technology 
Wright-Patterson  AFB,  Ohio 

.  attn:  asrmdf 

Aeronautical  Systems  Division 
Deputy  for  Systems  Management 
j  Support  Systems  Progrcuns  Office 
I  Wright-Patterson  AFB,  Ohio 
attn:  aSZT 

Armed  Services  Technical  Information 
Agency 

Document  Service  Center 
I  Arlington  Hall  Station 
^  Arlington  12,  Virginia  (20  copies) 

National  Aeronautics  and  Space 
Administration 
1512  H  Street  N.  W. 

I  Washington  25,  D.  C. 

I  attn:  Mr.  J.  Brewer,  Code  RA 

I. National  Aeronautics  and  Space 
Administration 
1512  H  Street  N.  W. 

Washington  25,  D.  C. 
attn:  Mr.  N.  F.  Rekos ,  Code  RAP 

National  Aeronautics  and  Space 
Administration 
Langley  Research  Center 
Langley  AFB,  Virginia 
attn:  Mr.  Donnelly 


National  Aeronautics  and  Space 
Administration 
Ames  Research  Center 
Moffett  Field,  California 
attn:  Mr.  W.  Cook,  40x80  Tunnel 

Office  of  Technical  Services 
Department  of  Commerce 
Washington  25,  D.  C. 

Library 

American  Institute  of  Aeronautics  and 
Astronautics 
2  East  64th  Street 

New  York  21,  New  York  (2  copies) 

Bell  Helicopter  Company 
P.  O.  Box  482 
Fort  Worth  1 ,  Texas 
attn:  Mr.  Robert  Lynn 

Cornell  Aeronautical  Laboratory,  Inc. 
4455  Genesee  Street 
Buffalo  21,  New  York 
attn:  Mr.  Frank  duWaldt 

Collins  Radio  Company 
Cedar  Rapids ,  Iowa 
attn:  Dr.  A.  Lippisch 

Georgia  Institute  of  Technology 
Guggenheim  School  of  Aeronautics 
Atlanta  13 ,  Georgia 
attn:  Mr.  D.  W.  Dutton 

Hiller  Aircraft  Corporation 
Advanced  Research  Division 
1350  Willow  Road 
Palo  Alto,  California 
attn:  Mr.  Raymond  Lockwood 

Massachusetts  Institute  of  Technology 
Aeronautical  Engineering  Department 
Cambridge  30,  Massachusetts 
ATTN:  Professor  R.  H.  Miller 

McDonnell  Aircraft  Corporation 
St.  Louis,  Missouri 
ATTN:  Dr.  Kurt  Hohenemer 

Mississippi  State  University 
Engineering  and  Industrial  Research 
Station 

State  College,  Mississippi 
attn:  Dr.  J.  J.  Cornish 


Naval  Postgraduate  School 
Aeronautical  Engineering  Department 
Monterey,  California 
ATTN:  Dr.  R.  Head 

Princeton  University 

Aeronautical  Engineering  Department 

James  Forrestal  Research  Center 

Princeton,  New  Jersey 

ATTN;  Professor  David  C.  Hazen 

Sikorsky  Aircraft  Division 
United  Aircraft  Corporation 
Stratford  1 ,  Connecticut 
ATTN:  Mr.  Philip  Michel 

Syracuse  University 
Mechanical  Engineering  Department 
Syracuse ,  New  York 
ATTN:  Dr.  S.  Eskinazi 

Therm  Advanced  Research  Division 
Therm,  Inc. 

Ithaca,  New  York 
ATTN;  Dr.  Ritter 

Vehicle  Research  Corporation 
1661  Lombardy  Road 
Pasadena,  California 
ATTN:  Dr.  Scott  Rethorst 

University  of  Virginia 
Aeronautical  Engineering  Department 
Charlottesville,  Virginia 
ATTN:  Dr.  G.  B.  Matthews 

Vertol  Division 

Boeing  Airplane  Company 

Woodland  Avenue 

Morton,  Pennsylvania 

ATTN:  Dr.  W.  Z.  Stepniewski 

Vidya  Division 
Itek  Corporation 
1450  Page  Mill  Road 
Stanford  Industrial  Park 
Palo  Alto,  California 
ATTN:  Dr.  J.  N.  Nielsen 


